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Abstract: For (3+l)-dimensional fermions, a net axial charge and external magnetic field 
can lead to a current parallel to the magnetic field. This is the chiral magnetic effect. We 
use gauge- gravity duality to study the chiral magnetic effect in large- N c , strongly-coupled 
N = 4 supersymmetric SU(N C ) Yang-Mills theory coupled to a number Nf -C N c of 
M = 2 hypermultiplets in the N c representation of SU(N C ), i.e. flavor fields. Specifically 
we introduce an external magnetic field and a time-dependent phase for the mass of the 
flavor fields, which is equivalent to an axial chemical potential for the flavor fermions, and 
we compute holographically the resulting chiral magnetic current. For massless flavors we 
find that the current takes the value determined by the axial anomaly. For massive flavors 
the current appears only in the presence of a condensate of pseudo-scalar mesons, and has 
a smaller value than for massless flavors, dropping to zero for sufficiently large mass or 
magnetic field. The axial symmetry in our system is part of the R-symmetry, and the 
states we study involve a net flow of axial charge to the adjoint sector from an external 
source coupled to the flavors. We compute the time rate of change of axial charge and of 
energy both in field theory and from holography, with perfect agreement. In contrast to 
previous holographic models of the chiral magnetic effect, in our system the vector current 
is conserved and gauge-invariant without any special counterterms. 
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1. Introduction 

Consider a (3+l)-dimensional system of free, massless Dirac fermions ip. The Lagrangian 
of such a system has two £7(1) symmetries, the vector one U(l)y, with conserved current 
Tp'y^tp, and the axial one U(1)a, with conserved current i/j-y^^tp. U(1)a is anomalous, 
and can be explicitly broken by a nonzero Dirac mass. If we introduce an axial chemical 
potential /is then we expect an imbalance in the number of left- and right-handed fermions. 
If we further introduce an external U(l)y magnetic field B then, assuming the fermions 
have positive charge, we expect their spins to align with B, and since they are massless their 
momenta will also align or anti-align depending on their chirality. Given the imbalance in 
chirality, we expect a net £7(1) y current parallel to B. This is the simplest example of the 
chiral magnetic effect (CME) [1-3]. 

For free fermions, the axial anomaly determines the size of the chiral magnetic current 
as follows. An axial chemical potential is equivalent to a background £7 (I) a gauge field with 
constant time component, = fj,$, or to a time-dependent phase ifi —> e^^il). Via the 
axial anomaly such a phase shift can be traded for a 0-term of the form a(t, x)F A F with 
a(t,x) = n^t and F the £7(l)y field strength. The spacetime-dependent source a(t,x) can 
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be regarded as a background, non-dynamical axion field. Writing F = dA and integrating 
by parts, we obtain an interaction of the form da A A A F. Varying the action with respect 
to A we find the chiral magnetic current, which is parallel to B and has magnitude 

J = 0F. (1.1) 

The quantity a = J/B is called the chiral magnetic conductivity. When the system has 
nontrivial time evolution a becomes a function of time, or in Fourier space a function 
of frequency a(oj). For free fermions, the DC limit <r(0) is fixed by the axial anomaly. 
Generically, interactions can modify a(oj), including <r(0) [4]. Notice that the CME occurs 
only when parity P and charge conjugation times parity CP symmetries are broken. 

A CME may occur in heavy-ion collisions such as those produced at the Relativistic 
Heavy-Ion Collider (RHIC) and the Large Hadron Collider (LHC) [2,5,6]. The dominant 
interaction in the early stages of collisions is the strong nuclear force, as described by 
Quantum Chromodynamics (QCD). Indeed, the quark-gluon plasma (QGP) created at 
RHIC appears to involve strongly-interacting degrees of freedom far from equilibrium: the 
plasma appears to thermalize quickly and have a short mean free path (both signs of strong 
interactions [7,8]) but also expands and cools rapidly until hadronization occurs. 

Two conditions must hold to produce a CME in a heavy-ion collision. First, the 
QCD vacuum apparently preserves P and CP, so some event-by-event violation of these is 
required. In a medium such as the QGP, one possible mechanism for such violations are 
fluctuations of the topological charge density. Second, the collision must be non-central, i.e. 
the nuclei must not perfectly overlap upon impact. In that case, the net charge combined 
with the net angular momentum can produce large magnetic fields, although these may 
die quickly as the QGP expands [2]. 

Assuming that P and CP are broken and a magnetic field is present, we know of two 
mechanisms to produce a CME in finite-temperature QCD. The first occurs for sufficiently 
large temperatures, where the QCD plasma is deconfined and chiral symmetry is restored. 
In that case we may invoke a naive picture of quarks as freely propagating fermions in a 
magnetic field, and apply the arguments above. 

A second, more subtle, mechanism, discussed for example in ref. [9], may occur at lower 
temperatures, when the QCD plasma is in a confined state with chiral symmetry broken. 
Here we expect a gas of hadrons rather than a QGP. The key observation is that an external 
electromagnetic field can convert a neutral pseudo-scalar meson, such as the ir°, rj, or rf, 
into a neutral vector meson, such as the p. More precisely, any effective action describing 
QCD and electromagnetic interactions will include for example a vertex of the form Bit® p. 
The vector meson so produced will be polarized in the direction of the magnetic field, and 
via interactions with charged mesons can induce a current parallel to B, thus producing a 
CME even in a confined phase. The same process may also occur in the late stages of QGP 
evolution, during hadronization when metastable domains with spontaneous breaking of P 
and CP could be formed [10]. 

To our knowledge, analysis of RHIC data appears to favor the presence of a CME in 
the QGP, although a better understanding of systematic errors and backgrounds is still 
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needed before a firm conclusion can be made [11—14]. The strong interactions and far- 
from-equilibrium evolution of the QGP in a heavy-ion collision make a clean theoretical 
prediction for a very difficult. Lattice simulations suggest that the CME occurs in thermal 
equilibrium [15-18]. Lattice simulations cannot yet reliably determine the time evolution 
of a, or equivalently a(oj), which is crucial for estimating the size of the chiral magnetic 
current in a heavy-ion collision. 

An alternative approach to the CME is the anti-de Sitter /Conformal Field Theory cor- 
respondence (AdS/CFT) [19-21], or more generally gauge-gravity duality. Gauge-gravity 
duality equates a strongly-coupled non-Abelian gauge theory with a weakly-coupled theory 
of gravity on some background spacetime, such that the field theory lives on the boundary 
of the spacetime, hence the duality is holographic. In particular, a black hole spacetime is 
dual to a thermal equilibrium state in which the center symmetry is spontaneously broken, 
such as the high-temperature, deconfined phase of a confining theory, where the tempera- 
ture of the field theory coincides with the Hawking temperature of the black hole [22]. 

Gauge-gravity duality has been most successful at describing out-of-equilibrium physics, 
especially near-equilibrium physics, i.e. hydrodynamics. Most importantly, all gauge the- 
ories with a gravity dual (in states with 5*0(2) rotational symmetry [23]) have the same, 
very small, ratio of shear viscosity rj to entropy density s, namely rj/s = 1/4tt [24], which 
is surprisingly close to the value estimated for the QGP at RHIC [25,26]. We take such 
universality, and indeed the universality of hydrodynamics in general, as encouragement 
to study the CME in many holographic systems, following refs. [27-33], including systems 
without confinement or chiral symmetry breaking in vacuum. 

A conserved U(l) current in the field theory is dual to a U(l) gauge field in the bulk 
and, roughly speaking, an anomaly for the current is dual to a (4+l)-dimensional Chern- 
Simons term for the bulk gauge field. The latter is thus typically a key ingredient in 
holographic descriptions of the CME [27-33]. More generally, holographic models dual to 
fluids with anomalous currents have been constructed for example in refs. [33-36]. These 
holographic studies are complementary to field theory studies of the effects of triangle 
anomalies on hydrodynamics [37-40], which themselves have been applied to study the 
CME in heavy-ion collisions [41-43]. 

One holographic model of QCD, the Sakai-Sugimoto model [44, 45] includes a bulk 
Chern-Simons term, although some confusion has arisen as to whether the CME occurs 
in this model at all. The problem in this model is that the vector current is anomalous 
under U(l)y x U(1)a transformations and therefore is not conserved in the presence of 
arbitrary external sources. Modifying the vector current such that it is conserved, which in 
the gravity dual requires adding certain boundary counterterms, causes the chiral magnetic 
current J to vanish [27]. To our knowledge, no consensus has emerged on whether a CME 
occurs in the Sakai-Sugimoto model. 1 

1 An alternative way to fix the normalization of the currents in the Sakai-Sugimoto model is to demand 
that the bulk action be invariant under gauge transformations that are non- vanishing at spatial infinity (in 
field theory directions) , which leads to different bulk counterterms [46] and produces a non- vanishing chiral 
magnetic current. For the sake of argument, here we are taking the phenomenological point of view that a 
U(1)a current in the presence of a U(l)v chemical potential should coincide with the weak-coupling result, 
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The authors of refs. [31,32] argued that for any bulk theory with a Chern-Simons term 
both a conserved vector current and nonzero J are possible, but at a price: the bulk gravity 
solution becomes non-regular. More precisely, in Euclidean signature the bulk gauge field 
will not vanish at the horizon and hence will not be a regular one-form [47]. In Lorentzian 
signature the gauge field solution will be regular only on the future horizon. One conclusion 
is that no reliable holographic description of the CME in thermal equilibrium (regular on 
the past and future horizons) exists for a conserved vector current. Effectively, a source 
for the gauge field must be introduced at the black hole horizon, the meaning of which is 
unclear from a field theory point of view. We were thus motivated to study other models 
where such issues could be avoided or at least clarified. 

We consider N = 4 supersymmetric SU (N c ) Yang-Mills (SYM) theory, in the 't Hooft 
large- N c limit and with large 't Hooft coupling, coupled to a number Nf <C N c of J\f = 2 
supersymmetric hypermultiplets in the fundamental representation of the gauge group, 
i.e. flavor fields. We introduce a complex mass m = \m\e 1 ^' for the flavor fields into 
the superpotential with a time-dependent phase <p = w ^ following refs. [48-51]. For the 
fermions that effectively introduces an axial chemical potential /X5 = ^oj. The theory also 
has a U(l)v symmetry that we will call baryon number. We introduce a baryon number 
magnetic field B and compute (holographically) the resulting chiral magnetic current. 

M = 4 SYM at large N c and large 't Hooft coupling is holographically dual to type 
IIB supergravity on AdS 5 x S 5 [19]. The N f < iV c hypermultiplets are dual to a number 
Nf -C N c of probe D7-branes extended along AdS§ x S 3 [52]. The phase of the mass 
corresponds to the position of the D7-branes in one of the transverse directions on the S 5 , 
hence to corresponds to the angular frequency of the D7-branes in that direction, and the 
axial charge density corresponds to the angular momentum of the D7-branes. The axial 
anomaly is realized holographically via the Wess-Zumino (WZ) coupling of D7-branes to 
the background Ramond-Ramond (RR) flux on the S 5 . The U(l)y current is dual to the 
U{1) gauge field on the worldvolume of the D7-branes, which thus encodes both the U(l)v 
magnetic field and the chiral magnetic current. 

In short, our system is a D7-brane in AdS§ x S 5 , rotating with angular frequency oj on 
the S 5 and with a worldvolume magnetic field B. At finite temperature we replace AdS§ 
with AdS-Schwarzschild. When |m| = the value of the chiral magnetic current agrees 
with the result from the calculation using the anomaly, eq. (|l.l| ). When \m\ is nonzero 
we find that a chiral magnetic current appears only when a certain [/(l)y-neutral pseudo- 
scalar operator has a nonzero expectation value, signaling the breaking of C times time 
reversal, CT. We interpret this as a neutral pseudo-scalar condensate being converted into 
a vector condensate by the magnetic field, in a manner somewhat similar to the CME in 
the low-temperature phase of QCD [9]. For nonzero |m| the value of the chiral magnetic 
current is less that that in eq. (1.1), and indeed both the current and the expectation value 



of the pseudo-scalar drop to zero for sufficiently large \m\ or B. 

Although we were motivated to find a model describing a CME in equilibrium with 
regular bulk solutions, in states with a nonzero |m| and a CME we can demonstrate that 



which occurs with the bulk counterterms of ref. [27] but not those of ref. [46]. 
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our system is out of equilibrium in two ways. First, we simply observe that the scalars in the 
N = 2 hypermultiplet have the same mass, with the same phase, as the fermions, so when 
\m\ is nonzero the Lagrangian has explicit time dependence and hence the system cannot 
be in equilibrium. The explicit time dependence disappears in the limit \m\ — > 0. Second, 
we observe that in our system the axial symmetry is part of the R-symmetry under which 
the adjoint fields of N = 4 SYM are also charged, hence axial charge in the flavor sector 
can "leak" into the adjoint sector, also taking energy with it. We compute, both in the field 
theory and from holography, the rate at which that occurs, with perfect agreement. We 
find that the rate is nonzero only when the pseudo-scalar has nonzero expectation value. 
Our solutions are stationary because we inject an equal amount of charge from an external 
source coupled to the flavor fields. The corresponding supergravity statement is that we 
pump angular momentum into the D7-branes which then flows into a bulk horizon. 

This paper is organized as follows. In section || we describe the main characteristics 
of the field theory with a flavor mass that has a time-dependent phase. In section ^ we 
describe the gravity dual and perform the holographic computation of the chiral magnetic 
current. In section || we compare with other holographic models and explain some of the 
issues related to current conservation, the N c dependence of anomalies, and to thermal 
equilibrium, and we compute the axial charge loss rate for our system. In section [5| we 
summarize and discuss our results. We collect some technical results in three appendices. 

2. The Theory in Question 

We study J\f = 4 SYM theory in the 't Hooft limit of iV c — > oo with Yang- Mills coupling 
squared g\ M — > keeping the 't Hooft coupling A = g\ M N c fixed, followed by the limit 
A> 1. The theory has an SO(6)r R-symmetry. The field content of M = 4 SYM theory 
is the gauge field, four Weyl fermions, and three complex scalars. The former are in the 4 
representation and the latter in the 6 representation of SO(6)r ~ SU(4:)r. We will also 
consider an N = 4 SYM plasma with equilibrium temperature T. 

We next introduce a number Nf of N = 2 supersymmetric hypermultiplets in the 
fundamental representation of SU(N C ), which in analogy with QCD we call flavor fields. In 
M = 1 notation the field content of the hypermultiplet is two chiral superfields of opposite 
chirality, Q and Q in the N c and N c representations, respectively. Each chiral superfield 
consists of a complex scalar and a Weyl fermion. We denote the scalars, the squarks, as q 
and q and combine the Weyl fermions into a Dirac fermion ip. The flavor fields' couplings 
break the SO(6)r symmetry down to 50(4) x U(1)ji, of which an SU(2)r x U (1)r subgroup 
is the Af = 2 R-symmetry. The U(1)r does not affect the squarks but acts as an axial 
symmetry for the quarks. Given that our theory has only this Abelian chiral symmetry, we 
will use "axial symmetry" and "chiral symmetry" interchangeably. As in QCD, the axial 
U(1)r symmetry is anomalous. The flavor fields also have a U{l)v symmetry that simply 
rotates Q and by the same phase. 

We will work in the probe limit, which consists of keeping Nf fixed when we take the 
't Hooft N c — > oo limit, and then working to leading order in the small parameter Nf/N c . 
Physically that corresponds to neglecting quantum effects due to the flavor fields, such as 
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the running of the coupling. For instance, when the 't Hooft coupling is small the probe 
limit consists of discarding diagrams with (s)quark loops. 

In the probe limit some part of the U (1)r anomaly survives, as we now explain. Three 
types of triangle diagram contribute to the anomaly, each with a U(1)r current at one 
vertex and two other currents at the other vertices. For example one diagram has the f7(l).R 
current and two gauge currents. We will denote that as the U(1)r x SU (N c ) x SU (N c ) = 
U(1)rSU(N c ) 2 anomaly, with similar notation for other anomalies. Both adjoint and flavor 
fields will appear in the loop, hence that diagram will have an order N% contribution and 
an order NfN c contribution. The next diagram is the U (1) R anomaly, which will similarly 
receive order N% and NfN c contributions. The third diagram is the U(1)rU(1) v anomaly. 
Only flavor fields carry the U(l)y charge, hence that diagram will be order NfN c , with no 
iV 2 contribution. In the probe limit we neglect the order NrN c contribution to the first 
two diagrams, since that is sub-leading. For the third diagram, however, the order NfN c 
term is leading, hence we retain it. 2 That anomaly will give rise to the CME in our system. 

N = 2 supersymmetry allows for a constant, complex mass m = \m\e 1 ^ for the flavor 
fields. A nonzero \m\ explicitly breaks U(1)r. We will introduce a mass with a time 
dependent phase: \m\e 1 ^ = \m\e tujt . Let us recall not only how <p = cot is equivalent to 
an axial chemical potential for the quarks, but also how <p = cot introduces explicit time 
dependence in the potential terms for q and q. Of the adjoint scalars, only one is charged 
under U(1)r. We denote this scalar as <I>. The flavor couplings in the j\f = 1 superpotential 
W are then 

W D Q$Q + \m\e i(t> QQ. (2.1) 

Integrating over superspace, we find the (normal) potential, from which we extract terms 
involving the squarks and terms involving the quarks. The terms involving q are [51] 

V q = q j \<S>\ 2 q - \m\e i4> q ] ^q - |m|e -< V$g + HV<7- ( 2 - 2 ) 

The potential includes identical terms for q. The quark contribution is simply 3 

= HV^V (2.3) 

If we now perform a chiral rotation 

Y> -> e-^ /2 V, (2.4) 

then the derivative in ^'s kinetic term will act on <f>, producing a new term that we may 
include in the potential, 

= |m|Y^ - M^ 7 5 ^. (2.5) 
If we introduce (j) = cot then clearly co is equivalent to twice the axial chemical potential, 

co = 2fi 5 . (2.6) 

2 On the supergravity side the first two anomalies would appear in the type IIB supergravity sector, 
while the third will be associated with a WZ term on a probe D7-brane. 

3 Our 7 5 is Hermitian and squares to the identity. 
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Crucially, when \m\ is nonzero the squark terms of the form in eq. (2.2) explicitly depend on 
t, so the Hamiltonian depends explicitly on time, energy is not conserved, and the system 
cannot be in equilibrium. 

For the CME in low-temperature QCD the central players are the light pseudo-scalar 
and vector fields, the 7r°, r/, rj and the p, respectively. Excitations of these fields, the 
mesons, produce poles in the corresponding retarded two-point functions. Our theory has 
operators analogous to these which will play a role in our realization of the CME, so let us 
describe them in detail. 

Our theory has gauge-invariant (s)quark bilinears, i.e. gauge-invariant operators built 
from two fields in the N c and iV c representations. The two-point functions of these operators 
exhibit poles which we will call mesons in analogy with QCD. Unlike mesons in QCD, these 
modes are not associated with chiral symmetry breaking or confinement, rather they are 
deeply bound states with masses on the order of \m\/^f\ [53]. When |m| is nonzero these 
are the lightest flavor degrees of freedom in our system. 

To determine the operators relevant for the CME, we treat \m\ and (j) as external 
sources. We denote the associated operators as O m and O^, respectively. For instance, 
varying minus the action with respect to |m|, we find a dimension three operator 

O m = --— = ^e^V - e^&q - e _i V$9 - e^tf&q - e~ irp q^q + 2\m\ ( q*q + tfq) . 
o\m\ V / 

(2.7) 

When (p is constant, O m is just the N = 2 supersymmetric completion of the standard 
quark mass operator. Notice that O m is charged under U(1)r, and hence may serve as an 
order parameter for chiral symmetry breaking when \m\ = 0. Notice also that if </> = u)t 
then O m depends explicitly on time. Varying minus the action with respect to 4> we fhrd a 
dimension four operator, 

= -— = \m\ iV^V^ + H q ] i (e" i0 $ - e^t) q + |m| q ] i (V^$ - q. 

(2.8) 

Notice that oc \m\, and again if (ft = ut then depends explicitly on time. 

The U(l)v baryon number and £7(1)^ currents will also be involved in the CME. We 
denote the conserved U(l)v current as J M , 

= - i (q^D^q - {D^q) ] q)-i(q {D^q) ] - (D^q) qA . (2.9) 



The contribution to the R-current from flavor fields is the same as half the axial current 
= ^V'7 /x 7 5 ' i / ; - 4 mentioned in the introduction, adjoint fields also contribute to the 
R-current, hence the axial current will not be conserved even in the absence of anomalies 
and when \m\ =0. We discuss the non-conservation of quark axial charge in detail in 



section 4.2 



4 We are identifying R-charge transformations with shifts <j) — > <j) + 8<t>, which for the quarks imply the 
U(1)a transformation il> — > e"^ 1 ^ 2 ip. With this convention the R-charge of the quarks is 1/2. 
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Since discrete spacetime symmetries play a central role in the CME, we will also present 
the transformation properties of various operators under C, P and T, when (f> = cot: 





C 


P 


T 


v q 


CO —7- — CO 


even 


CO — > —CO 




even 


CO — > —co 


even 




even co 


->■ —co + odd 


odd 


q \ i ( e »«t$t _ e -^«$) q 


co — > —co + odd 


even 


CO — > —CO 



where co — > —co means that a sign flip of co is the only change, and co — > —co + odd means a 
sign flip of co plus an overall sign flip are the only changes. The potential is not invariant 
under CPT, which is compatible with the breaking of Lorentz symmetry by the explicit 
time dependence. The only discrete spacetime symmetry under which the potential is 
invariant is CT. Notice that 0$ is CT odd, so an expectation value (0^) may serve as an 
order parameter for spontaneous CT breaking. 

Finally, let us explain the analogy between the CME in our system and in QCD. Our 
system has no dynamical electromagnetic U(l), so to obtain a CME we will introduce a 
non-dynamical external U(l)v magnetic field F xy = B. The phase <p will play the role of 
a non-dynamical external axion field a(t,x). The operator will play the role of a light, 
neutral pseudo-scalar, such as the 7r°. At zero temperature and finite mass, we can then 
think of an expectation value {O^) as a condensate of pseudo-scalar mesons. The U{l)v 
current will play the role of a vector meson field, like the p, so we can think of the chiral 
magnetic current (J z ) as a condensate of vector mesons. Our holographic calculations will 
show that (J z ) is nonzero only when (0^) is nonzero, except in the chirally symmetric case 
\m\ =0. Our interpretation is that the mechanism for the CME in our system is similar 
to that of low-temperature QCD: the magnetic field converts pseudo-scalar mesons into 
vector mesons polarized in the direction of B. Notice that away from the chiral limit the 
CME occurs in our system only when CT is spontaneously broken, in contrast to the free 
fermion case which required P and CP breaking. 



3. Chiral Magnetic Effect from Spinning Probe Branes 

We begin in type IIB string theory with a super symmetric intersection of N c D3-branes 
and Nf D7-branes: 





Xq 


Xl 


X2 


x 3 x 4 x 5 


Xq X 7 X 8 Xq 


D3 


X 


X 


X 


X 




D7 


X 


X 


X 


XXX 


X X 



Open strings with both ends on the D3-branes give rise at low energies to N = 4 
SU (N c ) SYM theory, while open strings with one end on the D3-branes and one on the D7- 
branes give rise to N = 2 hypermultiplets in the fundamental representation. The SO (6) 
isometry in the directions (X4, . . . , xg) corresponds to the SO(6)ji symmetry of N = 4 SYM 
theory. Clearly the D7-branes break that to SO{A) x U(1)r, corresponding to rotations in 
(X4, . . . 5X7) and (x8,xg) respectively. If we separate the D3- and D7-branes in the overall 
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= P (_^2 + d ^ + t_ ( dr 2 + r 2 ds 2 3 + dR 2 + (3 3) 



transverse directions, and xg, then the 3-7 and 7-3 strings acquire a finite length, giving 
the hypermultiplets a mass. The complex mass \m\e 11 ^ thus corresponds simply to the 
relative positions of the D3- and D7-branes in that plane, with |m| the separation distance 
and 4> the angle in the plane. The breaking of U{1)r by a nonzero \m\ appears simply as 
the breaking of rotational symmetry in the (xg, xg)-plane. A time-dependent phase <p = ojt 
corresponds to D7-branes spinning in the (xg, xg)-plane. 

We take the usual limits for the D3-branes, N c — > oo with <7 s iV c fixed, followed by 
taking <7 s iV c 3> 1, where g s is the string coupling and a' is the string length squared. We 
thus obtain type IIB supergravity in the near-horizon geometry of the D3-branes, AdS§ x S 5 
where each factor has radius of curvature L 4 /a' 2 = 47rg s N c » 1. The solution includes N c 
units of RR five- form flux on the S 5 . AdS/CFT equates this theory with the low-energy 
theory on the D3-branes, M = 4 SYM theory, with Yang-Mills coupling g\ M = ^g s an d 
't Hooft coupling A = gy M N c , so the theory is in the 't Hooft large- N c limit with A > 1. 

We will use an AdS§ x S 5 metric of the form 

ds 2 = -\gtt\ dt 2 + g xx dx 2 + g„ dr 2 + g S s ds 2 s3 + gRR dR 2 + g^ d<j) 2 , (3.2) 

p 2 

where p is the AdS^ radial coordinate with the boundary at p — > oo. The field theory lives 
in Minkowski space with coordinates (t,x). We have split the six directions transverse to 
the D3-branes into M 4 x M 2 where the latter R 2 represents the (x%, xg)-plane. We have 
written the metric of both the M 4 and the R 2 in spherical coordinates. The former has 
radial coordinate r with ds 2 s3 the metric of a unit S 3 , while the latter has radial coordinate 
R and circle coordinate 4>. Notice that p 2 = r 2 + R 2 . The self-dual RR five- form can be 
derived from a four-form potential 

C 4 = g 2 xx vol M s,i - g 2 ss dcf> A vol 5 3 . (3.4) 

Starting now we will use units in which L = 1. We can convert between string theory and 
supergravity quantities using a'~ 2 = A. 

The N = 4 SYM theory at finite temperature T is dual to supergravity in an AdS- 
Schwarzschild spacetime. In that case only \g tt \ and g xx change, becoming 

\9*\=P 2 y9$' 9xx = P 2 ^H(p). (3.5) 

with 

HP) = 1-^4, H(p) = l + -. (3.6) 

In these coordinates the horizon is always at p = 1, but the Hawking temperature, which 
we identify with the M = 4 SYM temperature, is T = 7/71", which we can vary by changing 
the parameter 7. We recover the T = limit by first rescaling p — > y/2p/j, and the same 
for r and R, and then taking 7—7-0. 

If we keep Nf finite as N c — > 00 we may treat the D7-branes as probes. The action de- 
scribing the D7-brane's dynamics, Sm, consists of two types of terms, a Dirac-Born-Infeld 
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(DBI) term and Wess-Zumino (WZ) terms. We will consider only the U(l) worldvolume 
theory of coincident D7-branes, so we will need only the Abelian D7-brane action, 



Sd7 = Sdbi + Swz, 



(3.7) 



Sdbi = -N f T D7 j d 8 ^ -det (g™ + (2ira>) F^ , 
Swz = +\n } T D7 (2to') 2 j P[C 4 ] A F A F, 



(3.9) 



(3.8) 



where To 7 




is the 



D7-brane tension, £ a are the worldvolume coordinates, g t 



is 



the induced metric on the brane, F a b is the U(l) worldvolume field strength, and PfCy is 
the pullback of the RR four-form to the D7-branes. 

Let us introduce some convenient notation. First, we will absorb a factor of (2ira') 
into the field strength (2ira r )F a i, = F a f,. Our D7-branes will be extended along AdS§ x S 3 
inside AdS$ x S 5 , that is, along the Minkowski coordinates (t, x), the radial direction r, 
and the S 3 C S 5 . In what follows we consider solutions for which the D7-brane Lagrangian 
will depend only on r, so we may trivially perform the integrations over the Minkowski and 
S 3 directions, producing factors of their respective volumes, Vjj 3 - 1 and 2ir 2 . We will absorb 
the factor of the infinite volume of Minkowski space into the action, Sdt/V^,i — > Sd7- 
From now on we will refer to this rescaled action as the D7-brane action. We will absorb 
the factor of the S 3 volume into an overall factor 



where in the second equality we converted to field theory quantities. 

We now need an appropriate ansatz for the worldvolume fields to describe a CME 
in the field theory. The two scalars on the D7-brane worldvolume are R and cf>. The 
former is dual to the operator O m while the latter is dual to 0<f>. More specifically, the 
asymptotic values of R and <p will be (proportional to) the modulus and phase of the 
complex mass, as is obvious from the initial D3/D7 intersection. We thus introduce R(r) 
and 4>(t,r) = ujt + ip(r), which produces a time-dependent phase for the mass and allows 
for nonzero |m|, (O m ) and (O^). 

We can motivate the r-dependence in (f)(t, r) from previous experience with probe 
branes in holographic spacetimes, following ref. [54]. Suppose we introduce only <j>{t) = oot 
and then perform a T-duality in the (j) direction. 5 The background solution changes, and 
the D7-brane becomes a D8-brane extended in <p with worldvolume gauge field, A^it) oc cot. 
The angular frequency oj becomes a constant electric field on the D8- brane, F t( f, a w. A 
probe brane in a gravitational potential well, such as AdS, and with a constant worldvolume 
electric field will generally have a tachyonic instability: the gravitational potential reduces 

5 Of course, T-duality is not a well-defined operation for an angular direction that can shrink to zero size, 
but here we are simply illustrating the similarities between spinning branes and branes with a worldvolume 
electric field. 



N = N f T D7 2?r 



XNfNc 



(3.10) 
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the effective tension of open strings, so at some point the constant electric field can rip 
strings apart. This tachyonic instability causes the Lorentzian action to become imaginary. 
Turning that around, an imaginary action signals a tachyon, since whenever a Lorentzian 
action S becomes a complex number, the weight factor e tS in a path integral will have 
either an exponentially growing or exponentially decaying mode. The cure for the D8- 
brane's instability is to introduce r-dependence [54], i.e. A<f,(t,r) oc wt + cp(r), producing a 
new constant of integration, associated with cp(r), that we can adjust to maintain reality 
of the action and hence avoid the instability. T-duality back to a D7-brane produces the 
4>(t, r) above, although now the physical interpretation of the instability is very different. 
Now the instability occurs because in a gravitational potential well the local speed of light 
decreases, while the probe brane rotates at a constant angular frequency, so at some point 
the probe brane may have linear velocity faster than the local speed of light. The D7-brane 
cures the problem by "twisting" in as a function of r [48-50] . 

For the CME we need a U(Y)y magnetic field and we expect a current {J z ). The 
U(l)v current J M is dual to the U(l) gauge field on the D7-branes, hence we include in our 
ansatz F xy = B and A z (r). In total, then, our ansatz includes 6 R(r), <fi(t,r) = ut + ip(r), 
F xy = B and A z (r). 

We can argue that the probe D7-brane action must depend only on derivatives of (ft and 
A z as follows. The background solution has an isometry in </>: the metric and four-form in 



eqs. (3.2) and (3.4) are invariant under constant shifts of 4>. Recalling that the scalars on the 
worldvolume of a D-brane are Goldstone bosons associated with the breaking of translation 
invariance in the transverse directions, and that Goldstone bosons can have only derivative 
interactions, we can conclude that the action Sd7 will involve only derivatives of (f>. For 
A z we argue simply that the action depends only on the field strength and not A z itself. 
Inserting our ansatz into the action, we find 



Sdbi = -M f drg^gH 2 J 1 + ?LJ (\g tt \ - g U 4> 2 ) (g rr + 9rr R' 2 + 9 XX A' 2 ) + \g a \g^' 2 
J y 9xx < 

Swz = -NBuj f drg 2 ss A' z . (3.11) 



where primes denote g- and dots denote As advertised, the action depends only on 
derivatives of <\> and A z and hence produces two "constants of motion," and s ^, 7 . 

We may thus solve for (j)' and A' z and obtain an action for R(r) only. We may do so in 
several ways. One is to solve for (j)' and A' z , derive i?(r)'s equation of motion from Sd7, and 
then plug the solutions for (j)' and A' z into that. Alternatively, we can plug the solutions 
directly into S07, perform a Legendre transform with respect to both <fi' and A' z , and then 
derive i?(r)'s equation of motion. We may also proceed in stages, for example by solving 
for and Legendre-transforming with respect to one only and then repeating the process for 
the second. The simplest approach turns out to be solving for eft' first and then for A' z . 

6 We work in A r — gauge. Recall also that our B includes a factor of (2ira!) oc A~ 1//2 . 
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The equation of motion for <fi' is 



SS D 7 Kr 3/2 3/2 /-, , B 2 \9tt\g<M><i>' 
~W = - M9 ^ dl V + gj 



a, 



■ ' 1 \gu\ - g^ 2 ) (grr + g RR R' 2 + g xx A' 2 ) + \g tt \g<i><t><t>' 2 

(3.12) 



where a is the first constant of motion. Solving for <fi' we get 

,2 (\gtt\- g<p<p4> 2 ) (g r r + gRR R' 2 + g xx A' 2 ) 



J2 



a 



~\gu\g H M 2 gl x gl s \g tt \g H (i + £)-o? 

Next we Legendre transform with respect to 4>' , 



(3.13) 



Sm = Sdbi + S W z = Sd7 - I drcf) ^TT- ( 3 - 14 ) 



Notice that Swz does not participate here, Swz = Swz, so we focus on Sdbi, 

Sdbi = S DB i - f drtf ^ S °f 1 (3.15) 



-N f dr g 3 s /2 g% 2 ^\g tt \ - g^ 2 ^ g rr + g RR R' 2 + g**A? J 1 + ^- 

j v gxx 



B 2 a 2 /M 2 



\gtt\guglxg S s 



The equation of motion for A' z is then 



5Sd7 5Sdbi . 8Swz a ,„ 

— + -YJT = ^ ( 3 - 16 ) 



SA' Z SA' Z 6A Z 



where /3 is the second constant of motion. The two terms in f3 are 



SSdbi 3/2 3/2 ^ ~2 9 XX A' Z B 2 a 2 /M 2 

= -Ng s ' s gj x \j \gtt\ - g H 4> 2 , = ==\l 1 + 



SA' Z V'**" ^ ^grr+gRRR' 2 + g xx A' 2 ^ g 2 x \gu\g H gl x gl s ' 

(3.17) 

^ = -MBo&s. (3.18) 

We now solve for A' 



A , (fi + NBujg 2 ss ) g xx yjg rr + gRRR' 2 

( M - (i + £ - c^u) - Z (f + AfW M 7' 

(3.19) 

Finally, we Legendre transform with respect to A' z , 

S D7 = S D7 - J drA'J-^ (3.20) 
= -M J dry / g rr + g RR R' 2 

I ~7~ ■ \ ( B 2 a 2 IN 2 \ f (3 \ 2 

x J AOs (m - (i + s - lgiiMxgls ) - *. ^ + ■ 
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We can derive R(r)'s equation of motion from this final form of the action, although is it 
cumbersome and unilluminating, so we will not present it. 

We can now explain how to extract field theory information from bulk solutions. The 
fields have the following near-boundary asymptotic expansions: 

rjf \ , c 2 , 1 2 lo g r , n f^gr\ 



a 



c 2 + \cquj 2 1 2 log r \ flogr 



Mr) = cz + o T n 4r- , (3.23) 



2 r 2 2 r 4 \ r 

where Co, C2, and c z are constants. 

In each case the leading term acts as a source for the dual operator, cq is the asymptotic 
separation between the original D3-branes and the D7-branes, so the magnitude of the mass 
is Co times the string tension, \m\ = j^j- The leading term in (f)(t,r), in our case cot, is 
the phase of the mass. c z is a source for J z , equivalent to the A z component of an external 
gauge field. In our case we may safely set c z = 0. 

The coefficients of the sub-leading terms determine the expectation values of the dual 
operators. The exact relations follow from the holographic dictionary, which equates the 
on-shell bulk action with minus the generating functional of the field theory. In appendix 
A we calculate the expectation values using holographic renormalization. For (O m ) we find 

(O m ) = (2vrQ')AA {-2c 2 - ^ 2 c - ^ 2 c In cjA . (3.24) 

Factors of the AdS radius, which we have set to one, make the argument of the logarithm 
dimensionless. 7 For the other operators we find 

(Otf,) = a, (J z ) = -(2ttc/)/3. (3.25) 

As shown in the last section, when \m\ is nonzero the operators O m and depend 
explicitly on time. In our solutions Co, C2, and a will be time-independent constants, 
however, so the expressions above are only consistent for nonzero \m\ if the state in which 
we evaluate (O m ) and {0<f,} has time dependence that cancels the time-dependence of the 
operators. Our configurations correspond then to a steady state. 

Notice that (O m ), {0$}, and (J z ) are all of different orders in the large- N c and large-A 
counting. (O m ) is of order (2ira')M oc ^/~\NfN c times factors of order one in the large-iV c 
and large-A counting, such as cq and C2- In contrast, (0</>) = a, where from eq. ( 3.12j) 



we see that a is of order J\f oc XNfN c times factors of order one, so {Oa,) is bigger than 
(O m ) by a factor of Recall, however, that is \m\ times a dimension three operator. 
Using \m\ = oc vAcq, we see that if {0$) scales as \NfN c then the expectation value 



At first glance, the (O m ) in eq. (3.24) appears to diverge in the flavor decoupling limit Co oc \m\ 00, 
which is counter-intuitive. Both the analytic argument in the appendix of ref. [49] and numerical calculations 
confirm that in fact (O m ) — >• in that limit. 
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of the dimension three operator must scale as (O t f ) )/\m\ oc \/XNfN c . The expectation 
values of the two dimension three operators, (O m ) and (0,jj}/|m|, thus have the same 
scaling. On the other hand we have (J z ) oc a '(3. From eq. ( |3T^ ) we see that in terms 
of large- A^ and large- A counting (3 oc NA' Z . Recall that we have absorbed a factor of 
(2-Ka 1 ) into A' z . If we extract that factor, then we find /3 oc {2na')N oc \^XNfN c and 
hence (J z ) = — (2-7ra')/3 oc NfN c . We thus find that (J z )'s normalization is independent 
of the coupling A. That is not surprising. (J z ) is our chiral magnetic current, whose 
normalization is fixed by the U(1)aU(1)y anomaly, and so is determined by the U(1)a and 
U(l)v charges, not the 't Hooft coupling. 

The massless limit |m| — > 0, or equivalently cq — > 0, is subtle. For one thing, in 
that limit the phase of the mass becomes ill-defined. Moreover, in that limit we see from 
cf)(t, r)'s asymptotic expansion in eq. ( |3.22| ) that a must also vanish, since otherwise the 
coefficients of the r-dependent terms in (p(t, r)'s expansion would diverge. The vanishing of 
a in that limit makes sense, since a = {O^) and we know from section |2|that (0^) oc |m|. 
As mentioned in section |2[ in our system (O^) is an order parameter for spontaneous CT 
breaking, so the vanishing of (£></,) as \m\ —> suggests that CT will always be restored in 
that limit. We must be cautious, however, since the expectation value of the dimension 
three operator (C,^) /| m| need not vanish as \m\ — > 0, so CT could still be broken. We will 
argue in section |3.2j that it actually is restored in the states we consider. 

Lastly, the R-charge density appears in the bulk as the angular momentum of the D7- 
brane. In what follows we will not compute the R-charge density or (O m ), rather our focus 
will be on (O^) and (J z )- We will only care whether the R-charge density is zero or not. 
To determine that we only need to know whether the embedding R(r) is nonzero or not: 
if R(r) = then the D7-brane has no angular momentum, while if R(r) is nonzero then 
the D7-brane has angular momentum. That is clear from the original D3/D7 construction, 
since if R(r) = then the D7-brane is not extended at all in the (x$, xg)-plane and so 
cannot have angular momentum. 

3.1 Solutions at Zero Temperature 



Our goal now is to solve -R(r)'s equation of motion, derived from eq. ( 3.2C| ), numerically in 



the pure AdS background dual to the zero-temperature vacuum of J\f = 4 SYM. 

Let us first quickly review what happens when B and oj are zero. Here we have no 
time-dependent phase for the flavor mass, and we expect no current, so we also set <fi and 
A z to zero. The induced metric on the D7-brane is then 8 

ds 2 D7 = p 2 (-dt 2 + dx 2 ) + ^(dr 2 (l + R' 2 ) + r 2 ds%), (3.26) 
where p 2 = r 2 + R 2 , and the action becomes 

S D7 = -M j drrVl + R' 2 - (3.27) 



8 Starting now we suppress the r dependence in R(r) for notational clarity, R(r) —¥ R, unless stated 
otherwise. 
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Inside the square root factor appearing in the action is a sum of squares, hence the action 
will be extremized only when R' = 0, or in other words when the solution is constant 
R = cq. These solutions describe flavor fields with an AT = 2 supersymmetry-preserving 
constant mass. N = 2 supersymmetry demands that (O m ) = 0, which is indeed the case 
for these solutions, which have C2 = and hence via eq. J3.24 ) (O m ) = 0. 



The D7-brane is always extended along r from the boundary r = oo to r = 0, however 
for these constant solutions the D7-brane does not fill all of AdS§. At the boundary the 
D7-brane wraps the maximum- volume equatorial S 3 inside the S 5 , but as it extends into 
AdS$, to smaller r, the S 3 shrinks and eventually collapses to zero size at the "North pole" 
of the S 5 , which occurs when r = 0. Recalling that the radial coordinate of AdS^ is not r 
but p = Vr 2 + R 2 , we see that at r = the D7-brane has only reached p = cq: from the 
perspective of an observer in AdS§ the D7-brane simply ends at that point. The trivial 
solution R = describes massless flavors. In that case the D7-brane fills all of AdS§. 

Notice that R = cq is a smooth solution because R'(0) = 0, that is, the slope of R is 
zero when the S 3 collapses at r = 0. If that does not occur then we see from eq. ( |3.26| ) that 
the D7-brane will have a conical singularity at r = 0. The regularity condition R'(0) = 
remains true when B and to are nonzero. In what follows we will find solutions for which 
R'(0) is nonzero and hence the D7-brane develops a conical singularity at r = 0. 

Now let us introduce a nonzero B, with <j) and A z still zero [55-58]. Roughly speaking 
a nonzero B "pushes" the D7-brane toward the boundary. More precisely, suppose we fix 
Co oc \m\ = 0. Here an infinite number of solutions appear, of which only one is the trivial 
solution R(r) = 0. The key question then is which solution has the smallest on-shell action 
and hence is physically preferred? A numerical analysis reveals that the trivial solution 
is not the preferred one: that honor is reserved for a nontrivial R(r) [55-58]. In fact, 
as we increase B the position where the D7-brane ends, p = R(0), increases. Physically, 
the nonzero B causes the D7-brane to "bend," and increasing B pushes the endpoint of 
the D7-brane closer to the boundary. Notice what that means in the field theory: the 
preferred solution, being nontrivial, necessarily has a nonzero C2, which from eq. ( 3.24| ) 



indicates a nonzero (O m ), hence chiral symmetry is spontaneously broken. The general 
lesson is that in our system a nonzero B promotes D7-brane bending, or in field theory 
language chiral symmetry breaking. The same remains true at nonzero cq , although in that 
case Co explicitly breaks chiral symmetry. 

Now let us return to nonzero B, (f>, and A z , and follow the arguments of refs. [48-51]. 
The induced D7-brane metric is now 



r 2 

ds 2 D7 = g% 7 dt 2 + Igf^dtdr + g^Jdr 2 + p 2 dx 2 + -^ds 2 s3 , (3.28) 



^ = P 2 ( "I + J^f) < ti 7 = 9?J = (3-29) 
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while the Legendre-transformed action in eq. ( 3.20 ) is 



S D7 = -M drr 2 Vl + R 2 (3.30) 



j 2 R 2 \ f B 2 a 2 /M 2 \ 1_ f l_ Bur 4 

(r 2 + R 2 ) 2 ) \ 1+ (r 2 + R 2 ) 2 ~ R 2 r e ) ~ r* \M + (r 2 + R 2 f 



The square root in the second line is of the form 9 a(r)b(r) — c(r) 2 where 

co 2 R 2 B 2 a 2 IN 2 1 (I , Bur 4 

a(r) " (r 2 + i? 2 ) 2 ' [r) ~ + {r 2 + R 2 ) 2 R 2 r 6 ' ° {r) ~ r 3 \M + (r 2 + R 2 ) 2 

(3.31) 

Notice that a(r) may change sign between r — > oo and r — > 0, but does not necessarily. 
More specifically, a{r) is always positive at r — > oo and may become negative as r — > 0, 
depending on the behavior of R(r). For the moment let us suppose that a(r) does change 
sign. We will denote the value of r where a(r) vanishes as ?v Upon taking a(r ! „) = and 
doing some algebra we find the equation for a semicircle, 

R(r*)-^Y+r 2 = ^ (3.32) 



2/ * 4 ' 

where the radius 10 is u/2 and the center is at (r*, i?(r*)) = (0, w/2). 

In fact, this semicircle is a horizon on the worldvolume of the D7-brane. If we change 
coordinates 

D7 

dt = dt + ^dr, (3.33) 
9tt 

then the induced metric becomes 

ds m = 9fi 7 di 2 + g^Jdr 2 + g xx dx 2 + gssds 2 s3 , (3.34) 

with 

D7 D7 D7 D7 (&? 7 ) 2 

i% = 5tt ) 9rr = Qrr WT~ " (3.35) 

We then have = —p 2 a(r) and hence a(r 4< ) = implies 9^ 7 {r*) = 0. To understand the 
appearance of this horizon on the D7-brane, consider a light ray moving in the cp direction, 
at fixed values of all other coordinates. The line element for a light ray is null, hence 
gttdt 2 + g^dcj) 2 = 0, which gives us the local speed of light in the 4> direction 



9ttl r+R (3.36) 



dt \ Qtheh R 



9 Notice that the Legendre-transformed action here has the same generic form as the Legendre-tranformed 
D7-brane action with worldvolume electric and magnetic fields used in ref. [59] for a holographic calculation 
of a Hall conductivity associated with the U(l)v symmetry. The similarity is not surprising, given the 
similarity between rotation and worldvolume electric fields due to T-duality, as explained above. Many of 
our arguments below are similar to those made in ref. [59]. 

10 Recall that we are using units in which the AdSs radius is L = 1. 
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Clearly when w is large enough to make a(r) < the D7-brane is moving faster than the 
local speed of light at that value of r, and a world volume horizon appears. Formally we 
can associate a temperature with the worldvolume horizon. We will discuss the meaning 
of this temperature in what follows, especially in appendix |C]. 

If a(r) changes sign but b(r) does not, then a(r)b(r) — c(r) 2 < for some r < r* 
and because of the square root Sd7 becomes imaginary, signaling a tachyonic instability 
as explained above. To avoid the instability we demand that b(r*) = also. Furthermore, 
as a(r) and b(r) approach zero, c(r) must approach zero more quickly, otherwise we again 
encounter an instability. We thus also impose c(r*) = 0. 

The condition a(r*) = fixes the worldvolume horizon while the conditions ^(r*) = 
and c(r*) = fix the two integration constants a and (3, or equivalently via eq. fl3.25| ) (O^) 
and (J 2 ). In other words, for given values of \m\, co, and B, regularity of the bulk solution 
determines unique values of the one-point functions (O m ), (0^), and (J 2 ), as is standard 
in AdS/CFT. Explicitly, we find 



/ B 2 Buir^ 

a = - mr ,)rl^ + - Jm - y d = -N mn) ; r , r (3.37) 

Using 0(7-*) = we can also express a and j3 in terms of R{r*) alone, 



a = -NR{n) 3 / 2 \R(r*)-u\ 3 /\lR(ny + ^r, f3 = -M-(R{n) - uof. (3.38) 



2 



The D7-brane does not always develop a worldvolume horizon. What happens when 
it does not? In that case a(r) > for all values of r. If a is nonzero then b(r) will change 
sign, rendering the action imaginary, so we demand a = 0. In addition if (3 is nonzero then 
again the action becomes imaginary because a(r)b{r) goes as 1/r 4 as r — > while c(r) 2 
goes like /3 2 /r 6 , so clearly \J a(r)b(r) — c(r) 2 will become imaginary at sufficiently small r. 
We thus also demand /3 = for these solutions. 

We thus have two classes of D7-brane embeddings, those with a worldvolume horizon 
and those without. The former have nonzero a and /3, or in the field theory nonzero (O^) 
and (J 2 ), hence these solutions describe a CME with CT spontaneously broken. The latter 
have a = and j3 = and hence no CME. 

An exceptional case is the trivial solution R(r) = and (p(t, r) = cot, corresponding to 
a chirally-symmetric state with \m\ = and (O m ) = 0. That solution has no worldvolume 



horizon, yet from eq. (3.38) we see that although a vanishes, f3 must be nonzero to maintain 
reality of the action. Recalling that B = (2Tra')B, where B is the value of the magnetic 
field in the field theory (see above eq. ( |3.10D ), we find for the trivial solution 

P = -{2na')NBoj. (3.39) 

Using eqs. ( |2.6D , ( |3.1C ), and ( p. 25 ) to translate to field theory quantities, we find 



(J z ) = ^F^sB, (3-40) 
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which is the value fixed by the anomaly, as expected in the chirally-symmetric case. We 
hasten to add three things. First, we will see that eq. ( |3.40| ) is unchanged at finite temper- 
ature, where the trivial solution remains a valid solution. Second, the trivial D7-brane has 
no angular momentum, so the corresponding field theory state has zero axial charge den- 
sity, despite having nonzero /is with zero mass gap, \m\ = 0. Third, because the anomaly's 
contribution to (J z ) does not contain much dynamical information, will will isolate the 
more interesting dynamical contributions by writing the (J z ) in eq. ( |3.25| ) as 

(J z ) = -(2irc/)P = {2na'fNBuo ( _ ) = ^^Boo ( _ ) . (3.41) 

V 1 V ' y ' \(2ira')NBojJ (2vr) 2 \(2^a')MBu: ) V ' 

The factor in parentheses in the final equality contains the non-trivial dynamical informa- 
tion. From eq. ( 3.37] ) we see that (3 oc (2ira')J\fBuj, so the factor in parentheses also does 



not depend explicitly on the magnetic field, although it will depend implicitly through the 
embedding. Notice that the current is always proportional to B, as we expect for the CME. 

To produce numerical solutions for the two classes of D7-brane embeddings, we must 
specify boundary conditions. The equation of motion for R(r) is a second-order non- 
linear ordinary differential equation, for which we need two boundary conditions on R(r). 
Solutions without worldvolume horizons are the simplest to produce. For these we set 
a = and (3 = 0, choose a value of R(0) greater than u/2 to avoid a worldvolume horizon, 
and then impose R'(0) = to guarantee regularity. Solutions with worldvolume horizons 
are trickier to obtain, 11 since the equation of motion itself depends on the values of a and /?, 
or equivalently on r* and R(r*), so we must choose these before we can solve the equation 



of motion. For these solutions we first choose a point on the semicircle in eq. (3.32), which 
fixes the values of a and (3 via eq. ( |3.3cf ) . We then obtain a condition on the first derivative 
at that point, R'(r*), from the equation of motion itself. We omit the explicit form, which 
is unilluminating. With these boundary conditions we can solve the equation of motion 
both inside the worldvolume horizon and outside. Notice that in these cases the value of 
-R'(O), which determines whether the D7-brane has a concial singularity at r = 0, is an 
output of the calculation. Figure |l] shows numerical solutions for R(r) for various values 
of \m\, u, and B. 

Our first observation is that for all solutions with a worldvolume horizon R'(0) is 
nonzero, so in our system at zero temperature all non-trivial solutions describing the CME 
have a conical singularity. In fact, for these solutions the on-shell action exhibits a diver- 
gence at r = 0, taking us outside of both the probe and supergravity limits, so strictly 
speaking we should not trust these solutions. Nevertheless, in section || we will argue that 
the singularities are physically sensible, being intimately related with the time rates of 
change of axial charge and energy in the field theory. 

Our second observation is that all solutions with a worldvolume horizon have nonzero 
I ml. That makes sense, since these solutions have nonzero a, and hence must have nonzero 



11 Solutions with nonzero B and worldvolume horizons but with A z (r) = were obtained in ref. [49]. 
These solutions are in fact unphysical, since an ansatz with A z (r) = is inconsistent: in ref. [49] the WZ 



term in eq. (3.11) was omitted, but the presence of that term necessitates the introduction of A z (r). All 



other solutions in ref. [49] besides these are consistent. 
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Figure 1: Numerical D7-brane embeddings R(r) for T — and lu = 1 for various values of B, in 
units of the AdS^ radius. The red semi-circle denotes the worldvolume horizon of eq. (3.32). (a.) 
B = 1. (b.) B = 5. (c.) B = 10. (d.) B = 20. The asymptotic value of R(r) as r -> oo (the far 
right in each plot) is the coefficient cq in eq. (3.21), which is proportional to the flavor mass \m\. 
The different classes of solutions, and their behavior as functions of B and \m\, are discussed in the 
accompanying text. 



\m\, as described above. Only solutions with a = can describe \m\ = 0. That class 
of solutions includes the trivial one R(r) = as well as non-trivial solutions without 
worldvolume horizons. 

From fig. [l] we can deduce the general behavior of solutions as we increase |m| or B, 
as follows. Suppose we fix u> and B, i.e. we choose one of figs, [j] (a.) through (d.), and 
then begin with some cq oc \m\ that is nonzero but much smaller than oj or \^B, such 
that the solution is very close to the trivial R(r) = solution. As we increase \m\, clearly 
the endpoint R(0) also increases. For sufficiently large \m\ the worldvolume horizon will 
disappear, at which point a and (5 vanish. Alternatively, suppose we fix ui and Co oc \m\, and 
then increase B. Now in fig. [l] we are choosing the value of a curve at the far right and then 
moving through the figures from (a.) to (d.). Again we see that for sufficiently large B the 
worldvolume horizon will disappear. 12 The corresponding field theory statements are that 
increasing \m\ or B eventually restores CT and extinguishes the CME, since eventually 

12 From fig. ^ (d.) we also see that for some D7-branes without a worldvolume horizon, R(r) passes 
through zero for sufficiently large B. Such behavior has been observed many times for D7-branes with 
worldvolume magnetic field: see for example refs. [55,58]. As argued in ref. [58], these solutions have a 
sensible interpretation in the field theory as a renormalization group flow, although in equilibrium they are 
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(a.) (b.) 

Figure 2: (a.) The value of (J z ), divided by the anomaly-determined value in eq. ( |3.40 ), as 
a function of the flavor mass divided by the axial chemical potential, |m|/^5. Here we set the 
magnetic field B = 1 and to = 2/is = 1 (in units of the AdS$ radius). The different curves 
correspond to different temperatures T — y/n, with green solid, yellow dotted, red dot-dashed, 
and blue dashed corresponding to 7 = 0,0.5, 1,2, respectively, (b.) The pseudo-scalar condensate 
(0$)/N versus |m|/jug, with B = 1 and u = 2/i 5 = 1, for the same temperatures as in (a.). 



{Off,) = and (J z ) = 0. The general lesson is that chiral symmetry breaking, whether 
explicit via \m\ or spontaneous via B, acts against the CME in our system. 

Our main result in this section is fig. ^. The green solid curve in fig ^ (a.) shows the 
exact behavior of (J z ), normalized to the value in eq. ( |3.40D , as we increase \m\/fi^, and 
and the green curve in fig. Q (b.) shows the same for (O^/ftf. At |m| = 0, (J z ) takes 
the value determined by the anomaly, while (O^) = 0. Increasing \m\/fi5, (J z ) decreases 
monotonically and eventually reaches zero, while {O^) increases, reaches a maximum, and 
then drops to zero. We omit the curves for (J z ) and (O^) versus V^B/^s, which are 
qualitatively similar to those in fig. ^. 

3.2 Solutions at Finite Temperature 

We now want to solve i?(r)'s equation of motion numerically in the AdS-Schwarzschild 
background, corresponding to an N = 4 SYM plasma at temperature T. 

Before doing so, let us briefly review what occurs when u = and B = 0, i.e. when the 
worldvolume gauge field and <j) vanish, summarizing refs. [60-65]. The main difference from 
the zero-temperature, pure AdS§ case is the presence of the AdS-Schwarzschild horizon, 
which divides D7-brane solutions into two categories. The first are similar to those in pure 
AdSs, namely D7-branes for which the S s shrinks and eventually collapses to zero size at 
some value of p outside the AdS-Schwarzschild horizon. The second category consists of 
D7-branes for which the S s shrinks but does not reach zero size by the time the D7-brane 
intersects the AdS-Schwzrzschild horizon. In the current context, solutions in the first 
category are called "Minkowski" embeddings while solutions in the second category are 
called "black hole" embeddings. In Euclidean signature, with compact time direction of 

not always the lowest-energy solutions. These D7-branes do not describe a CME and so are of less interest 
to us than D7-branes with worldvolume horizons. 
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period 1/T, the time circle collapses to zero size at the horizon. The two types of D7-brane 
solution thus have distinct topology: Minkowski embeddings have a collapsing three-cycle, 
the S 3 , while black hole embeddings have a collapsing one-cycle, the time circle. For 
Minkowski emeddings the condition to avoid a conical singularity when the S s collapses is 
R'(0) = 0, while the condition for black hole embeddings to avoid a singularity is that in 
the (r, R) plane the D7-brane must be perpendicular to the AdS-Schwarzschild horizon. 

When oj and B both vanish, the theory has one physically meaningful dimensionless 
parameter, T/\m\. Suppose we fix T/\m\ such that we have a Minkowski embedding and 
then increase T/\m\, say by holding \m\ fixed but increasing T. In the bulk the AdS- 
Schwarzschild horizon will grow and move toward the boundary, eventually encountering 
the D7-brane. The D7-brane solution then becomes a black hole embedding. Such a process 
involves a change in topology, so we have reason to expect that in general any observable 
associated with the flavor fields in the field theory will exhibit discontinuous behavior. 
Indeed, the bulk transition from Minkowski to black hole embedding appears in the field 
theory as a first-order phase transition [60-65]. 

Perhaps the most dramatic change in that transition occurs in the spectrum of D7- 
brane excitations, dual to the spectrum of mesons. For a Minkowski embedding the fluc- 
tuations of worldvolume fields are normal modes, i.e. standing waves trapped between the 
AdS§ boundary and the endpoint of the D7-brane. That translates into a field theory meson 
spectrum that is gapped and discrete [53]. For a black hole embedding the worldvolume 
fluctuations are quasi- normal modes, that is, the eigenfrequencies acquire an imaginary 
part. Physically, these fluctuations can leak energy into the AdS-Schwarzschild horizon 
and hence are damped. In the field theory the meson spectrum is gapless and continuous. 
The transition between the two is thus a kind of "meson melting" transition [66]. 

Introducing nonzero B, still keeping u = 0, qualitatively has the same effect as in 
the pure AdS§ case: increasing B pushes the D7-brane toward the boundary. If we start 
with a black hole embedding, for example, and keep T/|m| fixed while increasing B/\m\ 2 , 
eventually a transition occurs to a Minkowski embedding [55-58]. 

Now consider nonzero \m\, T, B, and uj. The induced D7-brane metric is then 

2 2 

ds 2 D7 = g% 7 dt 2 + 2g% 7 dtdr + g^Jdr 2 + p 2 ^-H{p)dx 2 + ^ds 2 s3 , (3.42) 

2 p 



where g^. and g^ r are the same as in eq. (|3.29| ) but 

^ 7 = P 2 f4^ + 7Z^LV /(p) = 1-4 H(p) = l + ' (3.43) 



v 2 H(p) (r 2 +i? 2 )V P 4 ' P 4 '' 

and we recall that p 2 = r 2 + R 2 , the AdS-Schwarzschild horizon is at p = 1, and the 
temperature is T = j/tt in our conventions. The location of the D7-brane's worldvolume 
horizon r* is now given by 

7 2 f(p*) 2 oj 2 R(n) 2 



2 H(p*) {rl + R{r*) 2 ) 2 



0, (3.44) 



or equivalently 



Y [ri + R{r*yy 
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Clearly the D7-brane worldvolume horizon is always outside of the AdS-Schwarzschild 
horizon, since /(p*) > /(p = 1) = 0. If w/7 oc — > then p* — )■ 1 and the two horizons 
coincide. At fixed mass, we may think of the [i^jT — > limit either as small at fixed T 
or large T at fixed /LI5. 

The physical arguments of the last subsection for the reality of the action are un- 
changed. Applying those arguments to fix a and (3 we find 



7 4 , / 4B 2 r 4 

(3.46) 

In AdS-Schwarzschild our D7-brane solutions fall into three categories. The first two 
are the straightforward generalizations of the categories of the last subsection: Minkowski 
embeddings without worldvolume horizons and Minkowski embeddings with worldvolume 
horizons. The new category consists of black hole embeddings, which necessarily have 
worldvolume horizons, as explained above. As in the last subsection, solutions without a 
worldvolume horizon describe field theory states with no CME and no spontaneous breaking 
of CT, while solutions with a worldvolume horizon, whether Minkowski or black hole, 
describe field theory states with a CME and spontaneous breaking of CT. 

The trivial solution R{r) = falls into the third category of embeddings, since it 



necessarily intersects the AdS-Schwarzschild horizon. From eq. ( 3.46| ) we see that for the 



trivial solution a = and /3 takes the value in eq. fl3~39p , so (J z ) again takes the value in 



eq. JOCD , hence we see that the value of (J z ) in the chirally-symmetric case, being fixed 
by the anomaly, is independent of temperature. 

For Minkowski embeddings the procedure to generate numerical solutions is the same 
as in the last subsection. In particular, for solutions with a worldvolume horizon we 
first choose a point on the worldvolume horizon and then use the equation of motion to 
determine the first derivative. We use the latter procedure for black hole embeddings too, 
since these necessarily have a worldvolume horizon. Notice that when we impose boundary 
conditions at the worldvolume horizon, the behavior of R(r) and its derivative at r = 
or at the AdS-Schwarzschild horizon, which determines whether the solution has a conical 
singularity, is an output of the calculation. 

Figure ^ shows numerical solutions for R(r) for various values of \m\, cj, B, and T. 
The results for Minkowski embeddings are similar to those of the last subsection. In 
particular, Minkowski embeddings with a worldvolume horizon have a nonzero R'(0) and 
hence a conical singularity. The black hole embeddings, however, do not have such a 
conical singularity: as figure || suggests, and numerical analysis confirms, in the (r, R) 
plane depicted the D7-brane "hits" the black hole horizon perpendicularly. 

Many of the conclusions from our T = analysis remain valid at finite temperature. 
All solutions with nonzero a have nonzero cq oc \m\. The worldvolume horizon eventually 
disappears as we increase \m\ or B: chiral symmetry breaking works against the CME in 
our system. 

Fig. |2| (a.) shows the chiral magnetic current (J z ), normalized to the anomaly- 
determined value in eq. ( |3.40| ), versus \m\/fj,5 for B = 1 and several values of T. At 
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B=l,(D=l,y=l 



(a.) 





(c.) 



(d.) 



Figure 3: Numerical D7-brane embeddings R(r) for ui = 1 and 7 = 1, corresponding to a temper- 
ature T = ~f/ir = 1/tt, for various values of 5, in units of the AdSs radius. The red quarter-circle 
represents the AdS-Schwarzschild horizon while the other red curve (the near-quarter-circle) de- 
notes the worldvolume horizon of eq. (|3.44 ), (a.) B = 0. (b.) 5 = 1. (c.) B = 10. (d.) B = 50. 
The different classes of solutions, and their behavior as functions of 5, |m|, and T, are discussed in 
the accompanying text. 



higher T the chiral magnetic current can persist to higher values of |m|///s before dropping 
to zero. Fig. ^ (b.) shows (O^) /N versus \m\/ ^ for B = 1 and the same values of T as in 
fig. H (a.). The qualitative behavior of the pseudo-scalar condensate is similar to the T = 
case, increasing, reaching a maximum, and then dropping to zero as we increase \m\/ix^. 
At higher T, the maximum of the condensate is larger, and the condensate also persists to 
higher values of \m\/ fi^. 

As mentioned above, seeing (O^) — > as \m\ — > is not enough to conclude that CT is 
restored in that limit. Oa, is \m\ times a dimension three operator, and the expectation value 
of that operator could remain finite as \m\ — > 0. For the states we consider, we can argue 
that the expectation value of the dimension three operator vanishes as \m\ — > as follows. 
In the limit \m\ — > we expect the solution to approach the constant one R(r) ~ cq oc |m|, 
and we expect the worldvolume horizon to approach the AdS-Schwarzschild horizon, so 
~ 1. Inserting these approximations into eq. ( 3.44| ) we find /(/?*) ^ \m\uj/^, and then 



from eq. ( |3.46 ) we find a oc \m\ 2 , so (O^) vanishes as \m\ 2 as \m\ — > 0, indicating that the 



dimension-three operator (C^)/|m| vanishes as \m\. We have confirmed that our numerical 
results for (O^) in fig. |2| (b.) behave as \m\ 2 as \m\ — > 0. 

For black hole embeddings we expect the spectrum of worldvolume excitations will 
be gapless and continuous, as in the uj = case. We expect the same for Minkowski 
embeddings with a worldvolume horizon: fluctuations of worldvolume fields will "see" the 
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worldvolume horizon as a genuine horizon, and hence we expect them to behave in a fash- 
ion similar to those of black hole embeddings. 13 More specifically, since we can associate 
a temperature with the worldvolume horizon, we expect the solutions for linearized fluc- 
tuations to translate into field theory two-point functions with a form characteristic of 
thermal diffusion at that temperature [51]. Moreover, in the bulk the worldvolume and 
AdS-Schwarzschild horizons generally will not coincide, so the worldvolume temperature 
will generally be different from the M = 4 SYM plasma temperature. This provides another 
hint that the system is not in equilibrium. 

Although we avoided an obvious instability by demanding reality of the D7-brane ac- 
tion, more subtle instabilities may exist in the spectrum of worldvolume fluctuations. In 
other words, some worldvolume fluctuations may still be tachyonic. Indeed, such instabili- 
ties have been found in very similar systems [69]. In appendix |B] we present a preliminary 
analysis of stability of linearized fluctuations. We find suggestive evidence that instabil- 
ities may occur for black hole embeddings with large enough magnetic field or chemical 
potential. We leave a complete analysis for the future. 



4. Comparing Holographic Models of the CME 

Our goal in this section is to compare our system to the previous holographic studies of 
the CME in refs. [27-33]. The first issue we discuss is holographic realizations of the 
U(1)a anomaly. For concreteness we focus on the Sakai-Sugimoto model [44, 45] because 
that involves probe flavor branes, just like our system, although much of our discussion 
applies for any theory in which the axial anomaly is realized holographically via a (4+1)- 
dimensional Chern-Simons term. The second issue we discuss is the regularity of bulk 
solutions describing the CME, which is intimately related to the question of whether the 
system is in equilibrium. Our states describing the CME are not in equilibrium, since axial 
charge and energy can leak to the adjoint sector. We compute the rates at which those 



occur in section 4.2. We end the section with proposals for new holographic systems that 



should exhibit a CME. 

4.1 Anomalies and Definitions of Axial and Vector Currents 

The Sakai-Sugimoto model is a holographic model of QCD that begins in type IIA string 
theory with the following non-supersymmetric intersection of N c D4-branes, Nf D8-branes, 
and Nf D8-branes 
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The low-energy theory living on the D4-branes is maximally-supersymmetric (4+1) di- 
mensional SU (N c ) YM theory. One spatial direction is compact, denoted (#4) above, 
with supersymmetry-breaking boundary conditions (fermions are anti-periodic in that di- 
rection). The holographic dual is type IIA string theory in the near-horizon D4-brane 

13 A gapless, continuous spectrum appears in the presence of a worldvolume electric field, which is closely 
related to rotation via T-duality as we argued above [67,68]. 
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geometry with compact X4 [22], which consists of an S 4 with N c units of four- form flux 
-F4 and a six-dimensional space with (3+l)-dimensional Poincare symmetry (at zero tem- 
perature), a radial/holographic direction, and the compact X4 direction whose size goes 
smoothly to zero at a finite radial position in the bulk. The radial and X4 directions thus 
form a "cigar" geometry. 

The supergravity approximation is only reliable in the 't Hooft large- N c limit with an 't 
Hooft coupling that is large at the X4 compactification scale. If the 't Hooft coupling were 
small at that scale, then the low-energy effective dynamics would be (3+l)-dimensional 
pure SU(N C ) YM theory: with supersymmetry broken, the SYM scalars acquire a mass 
due to loop effects while the fermions acquire a mass at tree level. When the 't Hooft 
coupling is large at the X4 compactification scale, the low-energy theory is SU(N C ) YM 
theory plus additional degrees of freedom, namely a tower of Kaluza-Klein modes associated 
with X4. 

The probe D8- and D8-branes localized in X4 introduce flavor fields [44,45]. At the 
(3+l)-dimensional intersection with D4 branes, the D8-branes introduce left-handed Weyl 
fermions while the D8-branes introduce right-handed Weyl fermions. Since we choose 
the same number of each, we can package them into massless Dirac fermions. The field 
theory then has a non-Abelian global flavor symmetry U(Nf)i J x U(Nf)n of which the 
U(1)l x U(1)r subgroup is anomalous. In the near-horizon D4-brane geometry the D8/D8- 
branes are extended along the S* 4 and the (3+l)-dimensional field theory directions, and 
describe a curve on the cigar. 

For the sake of comparison, we will consider also another holographic model of QCD, 
proposed in ref. [70], which uses the same compactified D4-branes but probe D6-branes 
instead of D8/D8-branes, 
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In fact, starting from our D3/D7 system, if we compactify X4 and then T-dualize in that 
direction we arrive at the D4/D6 intersection above. As in D3/D7, for D4/D6 the flavor 
fields are N = 2 hypermultiplets localized at the (3+l)-dimensional intersection, where the 
hypermultiplet scalars acquire a mass due to loop effects but chiral symmetry prevents the 
fermions from acquiring a mass. In the near-horizon geometry the D6-branes are extended 
along an S 2 C S 4 and the (3+l)-dimensional field theory directions, and also describe a 
curve along the cigar. Most of the features of the D3/D7 and D4/D6 systems relevant 
for the CME are the same, so we will treat these two systems on an equal footing in the 
following. For example, the D4/D6 system has only an Abelian chiral symmetry, just like 
the D3/D7 system. 

In the Sakai-Sugimoto and D4/D6 systems thermal equilibrium is realized holographi- 
cally by Wick-rotating to Euclidean time, compactifying the Euclidean time direction with 
period 1/T, and imposing anti-periodic boundary conditions for fermions in that direction. 
Two possible dual geometries are known, with the preferred one (dominating the bulk path 
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integral) depending on T. For T sufficiently small compared to the X4 radius, the bulk so- 
lution has a compact time direction but no horizon, and is dual to a confined state with 
the center symmetry preserved. At sufficiently large T a first-order phase transition occurs 
to a black brane geometry, dual to a deconfined state with spontaneously broken center 
symmetry [71]. 

The main difference between the Sakai-Sugimoto and D3/D7 or D4/D6 systems that is 
relevant for the CME is the bulk realization of the U (1) a symmetry. In the Sakai-Sugimoto 
model, axial charge is carried by the flavor fields only and not by the adjoint fields. In 
the bulk that becomes the statement that the dual U(1)a gauge fields are associated with 
the D8/D8-branes alone, and the axial anomaly is realized via Chern-Simons terms. More 
generally the U(1)l x U(1)r symmetries are dual to gauge fields Lm and Rm on the D8- 
and D8-branes, respectively. We define vector and axial fields as Vm = [Rm + £m)/2, 
Am = (Rm — Lm)/2- The U(1)a anomaly is then realized holographically via a WZ term 
in the D8-brane action 

Swz = ^(2ira') 3 N f T D8 [ P[C 3 ] A F A F A F, (4.1) 
a! JDS 

where Tjys is the D8-brane tension and P [C3] is the pullback of the RR 3-form giving rise 
to the four- form, F4 = dC%. Integrating by parts, using J S 4 F4 = N c , and adding the 
contribution from the D8-branes (corresponding to the right-handed Weyl fermions), we 
find 14 



qDS _ 
D W7, — ~ 



k J d 5 x e MNLP ^L M Fk L F^ Q + K Jd 5 x e MNLp QR M F* L FK Q , (4.2) 

where k = ^(2ira') 3 T£)8NfN c . We thus obtain (4+l)-dimensional Chern-Simons terms, 
which are gauge invariant up to boundary terms. In gauge/gravity duality that suggests 
the associated field theory gauge transformations are anomalous, which is indeed the case 
here: both U(1)l and U(\)r are anomalous. In terms of the axial and vector fields the 
Chern-Simons terms become 



S$z = l Jd 5 xe MNLp Q 



VmFV l F* q + l -A M FV L FV Q + 1 - Am F$ l f£ q 



(4.3) 



Both the vector and the axial combinations are anomalous, although the vector anomaly 
vanishes when F A = 0, i.e. when external axial field strengths vanish. 15 In the field theory, 
the first and second terms correspond to U(1)aU(1)y anomalies, while the the last term 
corresponds to a U(l)\ anomaly. In all cases the coefficient is of the same order, NfN, 
and survives the probe limit. 



Cl 



14 In this section, uppercase Latin letters M,N,... will denote (4+l)-dimensional bulk coordinates, in- 
cluding the holographic radial direction, while lowercase Greek letters /i, v, . . . will denote (3+l)-dimensional 
boundary field theory directions. 

15 Throughout this paper we only consider the effects of the anomaly on one-point functions of currents, 
i.e. on charge conservation equations. That is what we mean when we say an anomaly vanishes. The 
anomaly would still be visible in three-point functions of the current, even when external fields vanish, 
however. 
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In contrast, in the D3/D7 and D4/D6 systems the U(1)a is realized hole-graphically as 
rotations of the flavor branes in a transverse plane, which in the near-horizon limit becomes 
a transverse direction (ft inside the internal space, S 5 or S 4 . We have seen how that works 
for the D3/D7 system, so let us explain what happens in the D4/D6 system. The D6-brane 
action includes a WZ term 



where Tjjq is the D6-brane tension and F is the U(l) field strength on the D6-branes, dual 
to the U(l)v current. Given the similarity with the D3/D7 WZ term, we expect similar 
physics. The WZ term should be associated with the U(1)aU(1) v global anomaly. 16 The 
pullback of C3 will produce a factor of the derivative of <fi. In particular if we introduce 
dt<fi = uj and F xy = B, and if we demand a real on-shell action, then we expect a nontrivial 
A z and hence a CME in the dual field theory. 

Unlike the Sakai-Sugimoto model, the bulk gauge fields dual to U(1)a are not associ- 
ated with the flavor branes (D7- or D6-branes) alone. Instead the bulk U(1)a gauge field 
is a component of the metric that transforms as a gauge potential under diffeomorphisms 
of the form <p — > <p + Notice that a dimensional reduction to five dimensions can 

produce a Chern-Simons term for these gauge fields. Moreover notice that in order to 
introduce an axial field strength we must deform the geometry, i.e. we cannot use flavor 
brane worldvolume fields alone. 17 Bulk solutions with (from our point of view) an external 
axial magnetic field were constructed in refs. [73,74]. 

The shift <j> —>•</> + a, with constant a, corresponds to an R-symmetry transformation 
in the dual field theory: the U(1)a symmetry in the D3/D7 or D4/D6 systems is part 
of the R-symmetry, so in the field theory not only flavor fields but also adjoint fields are 
charged. This is a crucial difference with the Sakai-Sugimoto model, and has two major 
consequences. First, an axial charge in the flavor sector is not necessarily conserved, even 
in the absence of anomalies, since axial charge can leak into the adjoint sector. 18 Second, 
the adjoint fields can contribute to the axial anomaly. In the 't Hooft and probe limits 
the anomaly will have contributions of order N% and NfN c . We have seen that for both 
the D3/D7 and D4/D6 systems, the U(1)aU(1) v anomaly survives the probe limit, and is 
realized in the bulk via the probe flavor brane's WZ terms. 

Given the different realizations of chiral symmetry in these systems, chiral symmetry 
breaking will also be realized in different ways. In the Sakai-Sugimoto model in the low- 
temperature confining phase the D8- and D8-brane join, forming a U-shaped curve along the 
cigar, so that what appeared to be distinct U{Nf)i x U{Nf)n symmetries are broken to the 

16 Other WZ terms are associated with different anomalies, for instance the U(1)aSU(N c ) 2 axial anomaly 
was derived using a WZ term of the form f P[Gr] [72]. 

17 We cannot use the scalar (j> on the flavor brane worldvolume to introduce an axial field strength: as 
explained in section ^, the probe brane action will depend only on derivatives of <j>, i.e. on the exact 
one-form dip, which produces only zero field strengths d 2 (f> — 0. 

18 For the D4/D6 system in the confined phase, the dual field theory spectrum includes glueballs charged 
under the R-symmetry, so R-charge introduced in the mesonic sector alone need not be conserved. 




(4.4) 
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diagonal. That is interpreted in the field theory as spontaneous chiral symmetry breaking. 19 
In the D3/D7 and D4/D6 systems, on the other hand, chiral symmetry breaking means 
breaking rotational invariance in a transverse plane. We can do that either by separating 
the flavor branes from the color branes, corresponding to a flavor mass in the field theory 
and hence an explicit breaking, or by the flavor brane bending, as occurs for example when 
we introduce the magnetic field F xy = B (recall section 

Let us now compare in more technical detail the U (1)aU(1)y anomalies in both types of 
systems, with emphasis on the definition of gauge-invariant currents. We will schematically 
write the relevant parts of the probe brane actions as 

S m = S? 8 M (F L ,F R ) + K J d 5 xe MNL ^R M F§ L F p Q - k J d^xe MNL ^L M Fk L F^ 

(4.5) 

S D7 = S? 7 M {F V ) + K jd 5 x fl(r) e MNLP ^d M ^ L F^ Q . (4.6) 

The terms SBm and SyJ^ can be either the DBI action or a Maxwell's action obtained 
from an a' expansion to lowest nontrivial order, plus possible couplings to other fields. 
We actually don't care about the explicit form of these terms. We only care that they 
are gauge invariant. The factors of k for each system are straightforward to derive, being 
products of the brane tension, Nf and numerical factors. We have written Sd7 m the 
second line above, although the same form appears for the D6-brane action in the D4/D6 
system. We will use a radial coordinate r for both systems, with the holographic boundary 
at r — > oo. We have included a generic factor of the radial direction, O(r), in Sd7, since 
that may appear for general flavor brane embeddings in the D3/D7 and D4/D6 systems: 
for the D7-brane case see eq. ( |3.11 ). 



Crucially, notice that the actions in eqs. (|4.5| ) and ( |4.q ) are more general than the probe 
flavor systems we have been discussing. In fact the actions above are simply representatives 
of two categories of holographic systems describing the CME: those for which the U(l)v x 
U(1)a symmetry is realized via gauge fields (not necessarily arising from probe branes) 
with (4+l)-dimensional Chern-Simons terms and those for which the U(l)v symmetry is 
realized via a gauge field while the ^(1)^ symmetry is realized via rotation in <f). Most of 
what follows is relevant for any holographic system in one of these categories, although we 
will continue to use the language of probe branes. 

The equations of motion obtained from these actions are, in the Sakai-Sugimoto model 

- 9n^^ + ZKt MNLPQ F« L F« Q = 0, (4.7) 

- 9 N ^M- - 2,.e MNL ^Fk L F^ Q = 0, (4.8) 



19 A quark mass is possible in principle in the Sakai-Sugimoto model if the D8- and D8-branes are localized 
at the same point of the xa circle, in which case a tachyon field should be part of the low-energy description 
even close to the boundary. The tachyon action is not known in general, nevertheless this approach has 
been pursued in models based on the Sakai-Sugimoto model [75-77] and in an AdS/QCD context [78-80], 
where a scalar field typically describes chiral symmetry breaking [81,82]. 
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and in the D3/D7 or D4/D6 systems, 

ssBl 



-d N 



'YM 



Sd N V M 



+ Ane^^dN^ldL^Q = 0. 



(4.9) 



For the sake of argument, we will consider two different definitions of the field theory 
currents. The first comes from the equations for the radial components of the gauge 
fields R r , L r and V r , which take the form of conservation equations for currents in four 
dimensions, 



ZOD8 

0O u K r 



and 



xeD7 



o, 



0, 



0. 



(4.10) 
(4.11) 

(4.12) 



We will refer to the near-boundary r — > oo limit of the quantities in brackets above as the 
"conserved currents." In the r — > oo limit, the terms with explicit factors of the gauge fields 
R v and L u and derivatives in the field theory directions acquire their boundary values, so 
they become terms depending on the sources. 

The second definition of the field theory currents is the standard holographic one, 
from variation of the on-shell action. We will call these the "canonical currents." For the 
Sakai-Sugimoto model, varying R^ and separately we find 
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For the D3/D7 or D4/D6 systems we vary with respect to V M to find 
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ZOD7 



where we have used that Sym is gauge invariant, so for any of the gauge fields Am 

SSym &Sym 



Sd r A L 



(4.13) 
(4.14) 

(4.15) 
(4.16) 



We identify the canonical currents as the r — > oo limit of the quantities in brackets above. 
Notice that for the D3/D7 and D4/D6 systems the two definitions produce identical cur- 
rents. This is not true for the Sakai-Sugimoto model. 

In the Sakai-Sugimoto model the U(1)l and U(1)r conserved currents are not gauge- 
invariant, while the canonical currents are neither conserved nor gauge- invariant. The same 
is then true for the U(l)y and U(1)a currents, constructed as linear combinations of the 
U(1)l and U(1)r currents. We can add counterterms at the boundary that modify the 
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definitions of the currents, as discussed in ref. [27], which allow us to choose the vector 
current to be either conserved or gauge- invariant, but not both. This is the usual situation 
in field theory when the currents are anomalous. No matter how we choose to define the 
currents, the gauge- invariant part has an anomaly fixed by the equations of motion, so the 
vector current cannot be both conserved and gauge invariant. However, we can still make 
the vector current conserved and gauge-invariant under U(l)y transformations alone, i.e. 
we can eliminate the explicit appearance of Vm in the action in eq. (^^) by an integration 
by parts, 

S D8 = Ij d 5 xe MNL ^A M {^ L F V PQ + F$ L F$ Q ) . (4.17) 

The integration by parts produces boundary terms that can be canceled by a counterterm, 
or, turning things around, we can add counterterms that make the vector current conserved 
and gauge-invariant under U(l)v- Most important for the CME in the Sakai-Sugimoto 
model, with appropriate counterterms to make the U(l)v current conserved (though gauge- 
variant), the chiral magnetic current vanishes [27]. 

In the D3/D7 or D4/D6 systems, in the probe limit we have no explicit axial gauge 
field and so cannot introduce an external axial field strength. The vector current is then ob- 
viously gauge-invariant and conserved: the conserved current defined in eq. ( P^|) depends 
only on the field strength FX,. No special counterterms are required. 

The realization of the axial anomaly is somewhat more subtle than in the Sakai- 
Sugimoto model. The U{1)a symmetry is actually an R-symmetry, whose associated charge 
corresponds to angular momentum in (p. For our probe flavor branes that is intimately re- 
lated to the canonical momentum associated with the worldvolume scalar eft, which we 
denote n¥ , 



7T. 



M- SS D7 _m 



The second equation is just the <j) equation of motion. The probe flavor contribution to 
the expectation value of the R-symmetry current in the field theory directions, J^, is the 
integral of 7r^ over the worldvolume of the brane. 20 

W (4-19) 

Again, in the probe limit we have no external axial field strengths to contribute to an 
axial anomaly. The vector field strengths still contribute, however. If we distinguish two 
contributions to the angular momentum current from the YM action and the other from 
the WZ term, then the U(1)aU(1)y anomaly comes from the WZ contribution, 

(Jr)\wz = J dre r ^Q(r) F^ p . (4.20) 
Upon integrating by parts, we obtain 

(Jr)\ wz = -^ vap m Ki^Uoo + 4k / dre r ^d r n(r) V v F v ap . (4.21) 



Here and in section 



4.2 we ignore any potential divergences that may appear at the r — ¥ oo endpoint 



of the integral, which can be cancelled with counterterms that do not affect our main results. 



- 30 - 



Where we have dropped a term proportional to d u V r because it is a total derivative in 
the field theory directions. The second term on the right-hand-side is invariant under 
r-independent gauge transformations SV V = d v f{x^) since these also produce terms that 
are total derivatives in the field theory directions. Notice, however, the first term is not 
invariant under gauge transformations that are non-vanishing at the boundary. Adding 
the gauge-invariant YM contribution, we obtain the total probe flavor axial/R-symmetry 
current. The total current is not gauge invariant, due to the WZ contribution. Equivalently, 
we can define a gauge- invariant current that is not conserved, hence the axial/R-symmetry 
is anomalous. 

4.2 Thermal Equilibrium and the Holographic CME 

In the Sakai-Sugimoto model, with appropriate counterterms the U(l)y current can be 
made gauge-invariant and conserved, but the chiral magnetic current is then zero [27]. For 
any holographic model of the CME exploiting (4+l)-dimensional Chern-Simons terms, a 
nonzero chiral magnetic current can be achieved if the axial chemical potential is carefully 
defined, as argued in ref. [31]. In a bulk geometry with a horizon at some position r#, the 
gauge-invariant holographic definition of the field theory chemical potential is 



where in the second equality we chose A r = gauge. If At(rn) = then we may identify 
^5 = A t (oo). If Atirji) 7^ 0, then A t {co) is an external source, which is physically distinct 
from the chemical potential. Indeed, part of the argument of ref. [31] was that a nonzero 
^-t(rff) is necessary to obtain a nonzero chiral magnetic current. 

A serious issue arises with a nonzero A t (ru), however. If the system is in thermal 
equilibrium, then in the field theory we may perform an analytic continuation to Euclidean 
signature, where time becomes a compact coordinate with periodic boundary conditions 
for bosons and anti-periodic boundary conditions for fermions. On the gravity side we can 
analytically continue the metric to a Euclidean black hole, where the time coordinate is 
also compact and smoothly shrinks to zero size at th- The issue is that solutions with 
nonzero At(rji) are not regular in the Euclidean geometry: put briefly, the Killing vector 
4: vanishes at rjj, so the gauge field, which is a one-form, is regular only if At(rn) = 
(see for example ref. [47]). We may thus question the starting assumption, that the system 
is in thermal equilibrium. More generally, we may question whether gauge-gravity duality 
can describe a CME in a thermal equilibrium state at all. 

In fact, non-regular gauge field solutions appeared already in holographic studies of 
the hydrodynamics of charged fluids, including the effects of anomalies [34, 35]. 21 In these 
cases the bulk solutions represent small perturbations about static black hole solutions, 
including a perturbation of a bulk gauge field. The boundary conditions on the gauge fields 
are regular only on the future horizon: strictly speaking, the solutions do not represent 
a system in equilibrium, then perturbed, then relaxing back to equilibrium, but rather 

21 We thank Andreas Karch for explaining the following to us. 




(4.22) 
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solutions out of equilibrium in the infinite past that eventually settle into equilibrium in 
the infinite future. The key point is that non-regular boundary conditions on bulk gauge 
fields are a generic signal of out-of-equilibrium physics. 

Even when At(rjj) = 0, some holographic descriptions of the CME require non-regular 
solutions, if a time-dependent scalar is involved. Consider for example the realization 
of the CME in the soft- wall AdS/QCD model [29, 32], which is essentially AdS space 
with a certain profile for a dilaton-like scalar, representing the running of the field theory 
gauge coupling [82]. That model involves a (4+l)-dimensional Chern-Simons action, but 
also includes a term with the derivative of a scalar field, dual to a pseudo-scalar meson, 
analogous to the field <p on the flavor branes in the D3/D7 or D4/D6 systems (hence we 
use the same symbol, 0), 



where A is dual to the U(1)a current. As in the D3/D7 system, the CME appears when the 
time derivative of the boundary value of the scalar field is nonzero, dt4>(oo) ^ 0. Through 
a gauge transformation we can set A r = —d r (f), in which case the axial chemical potential 
becomes 



We can now set At(rjj) = and still obtain a nonzero chiral magnetic current, although 
regular Euclidean solutions should also have dtfyirn) = since 8m<P is a one-form. 

The upshot is that all of the holographic models of the CME that we know involve 
non-equilibrium physics. For the D3/D7 and D4/D6 systems we can also easily see that 
the CME occurs only out of equilibrium. The simplest observation is that the worldvolume 
scalar </> depends on time explicitly: dt<fi = oj. The equivalent field theory statement is that 
whenever \m\ is nonzero the action explicitly depends on time, as mentioned in section ^. 

More subtle non-equilibrium physics also occurs in the D3/D7 realization of the CME. 
Axial current conservation can be violated in three ways. The first way is explicitly via a 
nonzero \m\. The second way is due to anomalies. The third way is due to the fact that 
the axial symmetry is part of the R-symmetry, so that an axial charge density introduced 
in the flavor sector alone can leak into the adjoint sector. 

From the field theory point of view, axial charge can be lost to the adjoint sector in 
two ways, depending on whether mesons are melted or not. The equivalent bulk statement 
is that D7-brane angular momentum can be lost in two ways, depending on whether the 
D7-brane has a worldvolume horizon or not. 

If mesons are not melted, then the loss of R-charge occurs when mesons radiate R- 
charged glueballs. In the 't Hooft limit the interactions of color singlet mesons and glueballs 
are suppressed, so we expect the loss of R-charge to the adjoint sector not to be apparent 
in that limit. The dual statement is that a D7-brane with angular momentum but no 
worldvolume horizon can only lose angular momentum via radiation of closed strings, but 
the relevant interaction is proportional to the string coupling g s ~ 1/N C and so is suppressed 
in the classical supergravity limit, where g s — > 0. 




(4.23) 




(4.24) 
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If mesons are melted, then the spectrum in the flavor sector is no longer just delta- 
functions representing color-singlet mesons, but a continuum of modes whose interactions 
with the adjoint fields experience no large-iV c suppression. The dual statement is that 
a D7-brane with angular momentum and a worldvolume horizon (whether a Minkowski 
or black hole embedding) can transfer angular momentum across the horizon even in the 
super gravity limit. 

We will now confirm explicitly the above expectations by computing the time rate of 
change of axial/R-charge, both directly in the field theory and also from the D7-brane 
description, following refs. [51,83], with exact agreement between the two. Notice that 
when axial/R-charge leaks into the adjoint sector it should take energy with it, so we 
will also compute the time rate of change of the energy density, or more precisely the 
expectation value of the "tt" component of the stress-energy tensor (density), (T t t). 

We begin with the field theory calculations. For the rate of change of R-charge we just 
need to compute a Ward identity. Since an R-symmetry transformation is equivalent to a 
shift in the source, i.e. a shift in the phase of the flavor mass (j) — ¥ <f> -+- 5(p, we have, using 
the definition of in terms of a variation of minus the field theory action with respect to 
the phase of the flavor mass, 22 

^<J£> = (O+). (4.25) 

In our case nothing is changing in space, so we obtain for the time rate of change of the R- 
charge density dt(J R ) = {O^). Recall that oc |m|, so a necessary condition for the rate 
of change to be nonzero is for \m\ to be nonzero. In that case the potential and therefore 
the Hamiltonian are time-dependent, so the energy density will also not be conserved. 
Recalling the definitions of the potential terms V q and from section [2], the change in 
energy density is 

dt(T u ) = d t (V q + V^) = {Orfdth (4.26) 

where the second equality follows from the chain rule. In our case dt4> = w, and we just 
saw that (0$) = dt{J R ), so for our system 

d t (T tt ) =ood t (J R ). (4.27) 

The two rates of change are directly proportional. As expected, when R-charge leaks into 
the adjoint sector, it takes energy with it. 

We now turn to the bulk calculation of the same rates of change. The R-charge density 
(J R ) is given by the angular momentum of the D7-brane, 

poo poo rn 

<•*>-/„*"*-/„ *i£ (428) 

For concreteness we have written the lower endpoint of the r integration as r#, as appro- 
priate for a black hole embedding. For a Minkowski embedding the lower endpoint is r = 0. 
Taking dt of the charge density and using the (j) equation of motion in eq. (4.18), we find 

poo poo 

d t (J R ) = / dr fal = - dr d r 7r; = - ^ , (4.29) 



5„ 



22 The variation of the path integral is J d 4 x ( — i<9 fl <50?/>7 M 'y 5 i/; — 0<f,S(f>Y Integrating by parts and de- 
manding that the variation vanish for any 8(f), we find the relation in eq. (4.25). 
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where in the second equality we assumed homogeneity, so the derivative of iv¥ in any 
field theory spatial direction vanishes. Recalling from eq. ( 3.12| ) that -k t , = a, which in our 



system is independent of r, the R-charge density appears to be constant in time, (9t(Jjj) = 0. 
That is indeed true since the states we study are stationary. The reason why is nontrivial, 
however: the two terms in the final equality above cancel one another. The contribution 
from the lower endpoint represents the angular momentum that the D7-brane is losing, or 
equivalently the R-charge that the flavors are dissipating into the adjoint sector, while the 
contribution from the upper endpoint represents angular momentum that we are pumping 
into the system by hand via a boundary condition on the D7-brane, since we are forcing 
the D7-brane to rotate at the boundary, or equivalently in the field theory R-charge that 
we are pumping into the system from an external source. 23 Implicitly in our solutions we 
choose the latter precisely to cancel the former. The upshot is that the loss rate is the 
contribution from the lower endpoint, 

d t (J t R )\ loss = a=(0 4> ), (4.30) 

in perfect agreement with the field theory Ward identity. Only D7-branes with a world- 
volume horizon have nonzero a, hence in the field theory only states with melted mesons 
have a nonzero rate of change for R-charge (in the 't Hooft limit), in conformity with our 
field theory intuition. 

To compute the rate of change of the energy density, we need to compute the stress- 
energy tensor of the D7-brane. As explained in ref. [83], we can do that in two equivalent 
ways. The first way is directly, by variation of the D7-brane action with respect to the 
background metric. The second way is via a Noether procedure. Although we have used 
both methods, we will only present the latter, which is more efficient. Defining a Lagrangian 
via Sd7 = j drC, the stress-energy tensor density of the D7-brane, in the AdS§ directions, 
is 

The D7-brane stress-energy tensor is then the integral of 0^ over r. For values of M and 
N in field theory directions, we may equate the D7-brane stress-energy tensor's components 
with the flavor fields' contribution to the expectation value of the field theory stress-energy 
tensor [83], 

(Tt)= [ drSV (4.32) 

To compute the rate of change of energy density we will only need one component of the D7- 
brane stress-energy tensor density, Q r t = —uj^t = —oja. We proceed in a similar manner 
to the calculation of the R-charge rate of change. Using conservation of the D7-brane 
stress-energy tensor density, 8m @ A jv = 0> we find 

roo poo 

(dtT u ) = -(d t T\) = - / drd t e\ = / drd r @ r t = 9 r t |£ = -wa|~ . (4.33) 



23 From a field theory point of view, the external source does not have a clear, intuitive interpretation, 
in contrast to, say, dragging strings [84-86]. In those cases a stationary string solution represents a heavy 
quark being pushed through a plasma by a simple external source, a U(l)v electric field [84-86]. In our 
case we simply modify the potential by hand to produce the appropriate source. 
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The total energy is conserved, but the flux of energy at the boundary and the horizon is 
nonzero. The flux at the horizon corresponds to the rate of energy dissipation, 



d t (T tt }\ loss = u;a = u J (0 4> ), (4.34) 



where a is given by ( p6|) and is negative. More generally, the components @ r represent 



"external forces" acting on the probe [87]. The holographic calculation reproduces the 
relation between energy and charge loss rates, 24 

^m*)L s = ^(4>L s - (4-35) 

The dissipation rate is only of order a ~ XNjN c . That means that the flavor sector 
will only transfer an order N% amount of charge into the adjoint sector over a time of the 
order N c /\ oc l/gy M 3> 1. For times parametrically shorter than N c /X, we can ignore 
the dissipation rate and treat the background as a reservoir, in which case the stationary 
solution in the probe limit is a reliable approximation to the actual solution. This is 
similar to what occurs with constant electric fields on the D7-brane, where both energy 
and momentum (but not angular momentum) are dissipated in the bulk [83]. 

Recall that the above analysis is valid for Minkowski embedddings with worldvolume 
horizon (or for black hole embeddings in the zero-temperature limit) simply by taking 
th —7- 0. In section |3| we saw that such D7-branes have a conical singularity at r = 0. We 
can understand this singularity as a consequence of the angular momentum and energy 
flux along the brane. When the angular momentum and energy flowing along the brane 
reach the "bottom" at r = 0, they must be dumped into some source, or really a sink. 
For black hole embeddings that source is hidden behind the AdS-Schwarzschild horizon, 
and the part of the D7-brane outside of that horizon is non-singular. In the absence of the 
AdS-Schwarzschild horizon, and neglecting the backreaction of the D7-brane, the source is 
manifested as the "naked" conical singularity of the embedding. Something very similar 
occurs for the holographic dual of M = 4 SYM formulated on a spatial S* 3 with R-charge 
chemical potentials. Bulk solutions, known as "superstars" (dual to zero-temperature BPS 
states), exhibit naked singularities that have a sensible physical interpretation in terms of 
charged sources, namely giant gravitons [88]. 

4.3 Other Possible Holographic Models of the CME 

Finally, we will comment briefly on other possible holographic realizations of the CME, 
using probe flavor branes. The main objective is to devise systems where the probe brane 
WZ terms produce a coupling of the form d(j)f\F f\F ', although not necessarily from rotation 
in a transverse plane. We emphasize that these models are speculative, i.e. we have not 
actually constructed supergravity solutions for them. 



Consider D4-branes compactified on a spatial circle, as described in section 4.1. Instead 



of rotating probe branes, we can consider a black hole with angular momentum in the 



24 



For black hole embeddings the loss rates can also be extracted from (suitably regulated) divergences in 

5ui 



the angular momentum SSd7 and in the D7-brane stress-energy tensor at the AdS-Schwarzschild horizon 



as explained in ref. [83]. 
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internal space, or different fluxes in the background geometry. For instance, in the high- 
temperature black hole geometry we can introduce a background RR two-form F2 = dC±, 
where C\ = h(r)dt and h(r) approaches a constant, ho, as r — > 00. The meaning of this 
C\ in the field theory dual can be understood by introducing a probe D4-brane parallel to 
the color D4-branes. When the probe D4-brane is taken to infinity, the coupling to the C\ 
form becomes 

/ P[d] A F A F -> ho / A A F, (4.36) 

J DA J xyz 

where on the right-hand-side we have integrated over r. We can regard this as a Chern- 
Simons term in the effective three-dimensional theory that lives on the compactified D4- 
branes at high temperatures. We can thus associate C\ with a topological charge. Probe 
D6-branes couple to the C\ flux via a WZ term of the form J P[C\] A F A F A F. The 
D6-branes wrap an S 2 C S 4 , so we can see that a CME may be possible in the presence 
of magnetic flux 25 f s2 F^O. Similarly, probe D8-branes will have a WZ coupling of the 
form J P[Ci] AFAFAFAF, and hence may exhibit a CME if the worldvolume instanton 
number on the S 4 is nonzero, f s4 FAi 7 / 0. In the chirally-symmetric phase, the D8- and 
D8-branes are disconnected and cross the black hole horizon at different points [71]. In that 
case the instanton number on the D8-branes (tir) and D8-branes (til) can be different. 
The chiral magnetic currents for U{V)y and U(1)a will then be proportional to (with both 
parallel to the magnetic field), 

J v oz (n R - n L )B, J A oc (n R + n L )B. (4.37) 

The type IIB versions of the D4/D6 and D4/D8 intersections with nontrivial C\ are 
D3/D7 and D3/D9 intersections in the presence of a time-dependent axion Co- The relevant 
WZ term for the D7-brane is / P[dC } AAAFAFAF. A CME may be possible if we 
introduce non-zero flux on the internal three-sphere J sS A A F 7^ 0. The relevant D9-brane 
WZ coupling is J P[dCo] AAAFAFAFAF, and a CME may be possible with nonzero 
flux on the five-sphere f s5 A A F A F / 0. Notice that in these cases the boundary terms 
have been chosen in such a way that the probe brane actions are not gauge-invariant, so 
whether the U(l)v current is anomalous would need to be checked. 



5. Summary and Discussion 

We used AdS/CFT to study the CME in large- N c , strongly-coupled M = 4 SYM theory 
coupled to a number Nf <C N c of N = 2 supersymmetric flavor hypermultiplets. We intro- 
duced a time-dependent phase for the hypermultiplet mass, which for the hypermultiplet 
fermions is equivalent to an axial chemical potential, and we introduced an external, non- 
dynamical U(l)v magnetic field. When the magnitude of the hypermultiplet mass \m\ was 
zero, we found at both zero and finite temperature that the chiral magnetic current (J z ) 

25 A flux on the S 2 implies that the D6-branes are a magnetic source of RR four- form flux, which can be 
interpreted as a bound state of D6-branes with D4-branes in the directions transverse to the sphere. Similar 
statements apply for the D8-, D7- and D9-branes we discuss next: flux on the internal space represents a 
bound state with D4- or D3-branes, respectively. 
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coincided with the weak-coupling result in the chir ally-symmetric limit, eq. ( |3,40D . When 
| to | was nonzero we found that (J z ) had a smaller value than eq. ( 3.40| ), and also that the 
J7(l)y-invariant and CT-odd pseudo-scalar operator acquired a nonzero expectation 
value. Indeed, for sufficiently large \m\ or B, compared to the axial chemical potential 
or the temperature, both (J z ) and (0^) dropped to zero (recall fig. Q). In these cases we 
interpret the appearance of a chiral magnetic current as the conversion of the pseudo-scalar 
condensate to a vector condensate via the magnetic field. 

We compared several holographic models of the CME, and in particular highlighted 
the fact that all of them describe a CME only in non-equilibrium states (in contrast to 
lattice QCD analyses [15-18]). Our holographic system also describes the CME only in 
non-equilibrium states. Whenever \m\ is nonzero, the scalars in the hypermultiplet have 
masses with time-dependent phases, hence the Hamiltonian has explicit time dependence 
and energy is not conserved. Moreover, in our system the axial symmetry is part of the 
R-symmetry, so axial charge in the flavor sector can leak into the adjoint sector, also taking 
energy with it. We computed the associated loss rates, which we found to be proportional to 
{0<h). When the CME occurs in our system at nonzero \m\, the flavor fields are losing axial 
charge and energy to the adjoint sector. In the probe limit these loss rates are negligible, 
however, and our states describing a CME were in fact stationary. 26 

The supergravity description of the above was a number Nf of probe D7-branes ex- 
tended along AdSz, x S 3 inside AdS$ x S 5 , rotating on the S 5 and with worldvolume gauge 
fields that encode the magnetic field and chiral magnetic current. AdS space is effectively 
a gravitational potential well in which the local speed of light decreases as we move away 
from the boundary. A D7-brane rotating sufficiently quickly may at some point be ro- 
tating faster than the local speed of light and hence may develop a worldvolume horizon. 
We saw that indeed D7-brane solutions thus split into two categories, those that rotate 
quickly enough to develop a worldvolume horizon and those that don't. For the former 
the D7-brane becomes imaginary, signaling the presence of a tachyon, unless we introduce 
certain worldvolume fields and adjust their integration constants to maintain reality of the 
action. Via the holographic dictionary these integration constants were precisely the values 
of {J z ) and {O^). In the bulk the loss of axial charge and energy appear as the flow of 
angular momentum and energy across the worldvolume horizon. We found numerically 
that the angular momentum and energy flux produces a conical singularity in Minkowski 
embeddings with worldvolume horizon at the point where the S 3 collapses to zero size. 

Although we focused on D7-branes, the CME can be realized in many similar flavor 
brane systems. The basic ingredients are a holographic spacetime with probe flavor D- 
branes satisfying two conditions: they describe (3+l)-dimensional flavor fields, and they 
have at least two transverse directions in which to rotate. An axial chemical potential, 
implemented as a time-dependent fermion mass, will be realized via rotation in a transverse 
plane, and the axial anomaly will be realized via a WZ coupling to RR flux in the internal 



26 The loss rates are of order (0^) oc \N C , and so can be neglected for times shorter than N c /X. Taking 
into account the change in angular momentum and energy, i.e. computing the back-reaction of the D7- 
branes, would probably lead to an expanding horizon [51]. Gravity solutions exhibiting precisely that 
behavior have been constructed for external electric fields in ref. [891. 
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space. A model relevant for applications to QCD is that of ref. [70], with flavor D6-branes 
in the near- horizon geometry of D4-branes. 

From a phenomenological point of view, models of the D3/D7 or D4/D6 type have ad- 
vantages and disadvantages when compared to the Sakai-Sugimoto or AdS/QCD models. 
Consider first the disadvantages. D3/D7-type models are typically less similar to large- 
N c QCD, for example, non-Abelian chiral symmetries will generically be explicitly broken 
by (super)potential terms. Given that the general objective of holography is to uncover 
universal physics, this disadvantage may not be fatal. Beyond that we suspect that the 
problems with D7-branes will be generic: U(1)a charge may leak into the adjoint sector, 
Minkowski embeddings describing a CME will likely be singular, and many solutions de- 
scribing the CME may in fact be perturbatively unstable. On the other hand, the fact 
that Sakai-Sugimoto or AdS/QCD models require a source at the horizon to describe the 
CME suggests that the same (or similar) problems may appear in those models as well, 
once the meaning and effects of the source are clarified. Perhaps the principal advantage 
of D3/D7-type models is that the U(l)v current is conserved and gauge invariant under 
U(l)v transformations by construction, so among other things comparison with weak cou- 
pling calculations is more straightforward. Another advantage is that quark masses are 
easy to introduce and the effects of chiral symmetry breaking are easy to study. 

Although our focus was on the CME, our spinning D7-branes without worldvolume 
horizon may have useful applications as well. When the solution describes massless flavors, 
the field theory is in a state with a finite charge density of fermions with U(1)a sponta- 
neously broken, i.e. superfluid states [49]. These solutions have no loss of axial charge or 
energy and no obvious instabilities, at least to leading order in the 1/N C expansion, and so 
deserve further study as models for strongly-coupled, many-body fermion physics. 27 

An important task for the future is a complete analysis of linearized fluctuations of 
worldvolume fields, to determine whether our solutions are stable. In particular, in our 
analysis we assumed homogeneity of the ground state, but in QCD with U(1)a or U(l)y 
chemical potentials and a strong magnetic field, in a state with chiral symmetry broken, the 
ground state may be the inhomogeneous "chiral magnetic spiral" [90]. Such a phase was 
indeed detected in the Sakai-Sugimoto model via analysis of linearized fluctuations [91], 
and a similar analysis should be done for the D3/D7 model. 
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A. Holographic Renormalization 



In this appendix we compute holographically the one-point functions of the operators O m , 
and J z defined in section ||. 

The AdS / CFT dictionary equates minus the on-shell bulk action with the CFT gener- 
ating functional. The on-shell action is generically divergent, however, due to the integra- 
tion over the infinite volume of AdS space. In the field theory these are UV divergences 
which require renormalization: we introduce a regulator, add counterterms, take varia- 
tional derivatives of the generating functional to obtain regulated correlators, and then 
remove the regulator to obtain finite, renormalized correlators. The analogous procedure 
in AdS /CFT is called holographic renormalization (see for example refs. [92,93]). Holo- 
graphic renormalization for probe D-brane worldvolume fields has been studied in various 
places. We will follow refs. [49,94,95]. 

We begin by writing the metric of AdS^ xS 5 (in L = 1 units) as 



ds 2 



1 



~2 (du 2 - dt 2 + dx 2 ) + d6 2 + sin 2 6d(j) 2 + cos 2 Odsgs, 



(A.l) 



where u is the radial coordinate with the boundary at u = 0. These coordinates are related 

(A.2) 



to those of eq. (3^) as 



1 



1 



u 



9 = arctan 



R 



p 2 r 2 + R 2 ' r 
The ansatz for the worldvolume fields described in section || is now 

6(u), A z (u), A y (x) = Bx, 4>(t, u) = uit + ip(u). 



(A.3) 



The metric and ansatz are identical in form to those of section ||, so formally we may 
use all of the equations of that section, with the replacements r — > u and R(r) — > 8(u). 
In particular, we may derive #(ii)'s equation of motion from eq. (3.20), with r — )■ u and 
R(r) -> 6(u), 



N / du\J g uu + 



x \ 9xx9ss (\9tt\ 



1 + 



£?2 

9xx 



a- 



N 2 \gtt\g<f><f>gL9ss 



1L 



+ Bu} 9ss 
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From the equation of motion we find 0(u)'s asymptotic expansion, 



6{u) = Cq u + &2 u 3 



u 3 log u + O (u 5 log u) . 



(A.4) 



A straightforward exercise using eq. ( |A.2| ) shows that the coefficient of the leading term, co, 
is identical to that of R(r)'s expansion in eq. fl3,21| ), hence we use the same symbol. The 
coefficients of the sub- leading terms in eqs. (EOlD and flAgj) are related as — C2 + ^Cq. 
For the rest of this appendix we use only cq and C2 unless stated otherwise. 

From eqs. ( p. 13 ) and ( 3.19| ), with r — > u and R(r) — > 9(u), we obtain the asymptotic 
expansions of the other fields 



A z (u) 



a 



1 



W, -)=^ + ^{^ + ^--\ C ^-^J^ + - 4 ^uHo gU ) + 0^ logn) . 



1 u 2 

=0 



lc 2 



We now plug 0(u), 4>(t,u) and A z (u) into Sd7 to find the divergences of the on-shell 
action. To regulate the divergences we integrate only to u = e, producing the regulated 
action 



qreg 



J\f I du 
1 1 



4e z 



1 



+ 



+ 



B 2 - 2oj 2 4 
2u 



+ ... 



r 2 1 

ro _ 

2 e 2 



2 cgloge + - J B 2 loge + 0(e 2 log e ) 



+ .... 



(A.5) 
(A.6) 



where we have presented only the terms that diverge as e — > 0, with the . . . representing 
all other terms. 

We next add counterterms at the u = e surface to remove the e — > divergences. The 
counterterms we need are those written in refs. [94,95] for the field 6(u), but modified 
to account for (p(t,u) as follows. Roughly speaking, 9{u) is dual to the magnitude of the 
mass operator, while <f> is dual to its phase. The precise operators are written in section |2[ 
When both scalars are active, the natural thing to do is thus to package these two fields 
into a single complex scalar field 0(i, u) = 9(u)e' l ^ t ' u \ and rewrite the counterterms of 
refs. [94,95] in terms of Q(t,u): 



L\ = +J\f ■ 



La = +M- 



l 



-7|0| 



-7 iog|e|e*n 7 e, l 5 = +m 



1 r 



7|9| 4 , 



12 

~e*n 7 e, 



(A.7) 
(A.8) 



U = ~ M ^ v 7 ^ log eF^F^, 



(A.9) 



its determinant, and 



where = e 2 "q liu is the induced metric at u = e, with 7 = — 
□~ is the scalar Laplacian associated with 7^^. Inserting the asymptotic expansions for 
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the fields, we find 



Ll = +M \^ (A - 10) 

L 2 = +AT -c c 2 + ^lo 2 4 log + 0(e 2 loge), (A.ll) 

L 3 = +M^c 4 + O(e 2 loge), (A.12) 

L A = +N (+\u 2 Co log e + ^u 2 c 2 log c ^ + O (e 2 log e) , (A.13) 

L 5 = +M^u; 2 c 2 + 0(e 2 loge), (A. 14) 

L e = +AA (-^B 2 \oge \ . (A.15) 



Comparing the above equations for the counterterms with the regulated action in eq. ( |A.6| ), 
we see that all divergences will cancel. 

Notice that L3 and L5 are finite in the e — > limit. We could also introduce a third 
finite counterterm, identical to Lq but without the factor of log e. Fixing the coefficients 
of these finite counterterms corresponds to a choice of renormalization scheme in the field 
theory. The coefficient of the |G| 4 counterterm is fixed by demanding that the one-point 
function (O m ), which we will compute shortly, vanishes for the solution that represents a 
constant J\f = 2-supersymmetric flavor mass at zero temperature with B = and uj = [94] . 
The coefficient of the 0*D 7 counterterm is fixed by demanding that (O m ) vanish in the 
large-mass limit \m\ cx cq —> 00, as demonstrated in the appendix of ref. [49]. For the 
third finite counterterm, proportional to B 2 , we follow the convention ref. [96] and set the 
coefficient to zero. 

Now let us compute renormalized one-point functions. To do so we define a "subtracted 
action" as 

S5* = s$ + X>. (A.16) 

1=1 

The finite, renormalized on-shell action is then 

Sffi = )imS$. (A.17) 

The AdS/CFT dictionary then equates —Spy with the generating functional of the field 
theory. For a (pseudo) scalar field <&(v) dual to an operator O of dimension A, the renor- 
malized one-point function is 

(P)=-lim^ ' fjl . (A.18) 

Our scalar 0(u) is dual to the dimension-three mass operator of the flavor fields, O m . More 
precisely, O m is given by taking the variational derivative of the generating functional with 
respect to \m\. Recalling from the discussion below eq. ( |3.23| ) that \m\ = t^Jj we find 

1 1 c J sub 

(O m ) = -(2W) lim _-__2Z.. (A.19) 
e^o e A y/^i 59(e) 
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The calculation of (O m ) is straightforward, so we present only the result, 
(O m ) = (2na')M {-2~c 2 + ^4 - ^ 2 c - ^ 2 c lnc^ . 



(A.20) 



Upon using C2 = C2 + gCg we obtain eq. ( |3.24| ), which is identical to the result for (O m ) in 
ref. [49]. The pseudoscalar <fi(t,u) = uot + (p(u) is dual to the dimension-four operator 
written in eq. (2.8), so for (O^) we have 



1 1 < ? SMb 



Notice that 



e->o e 4 y /:z 7 5<p(e) ' 
= = 1 so we only need to compute the variation, 

r x osub r 

SSj# = / du —^-d u S(p = a / dud u 5ip = -a5tp(u 

J 6 M uV Je 

SSf 



(A.21) 



(A.22) 



where we used the fact that on-shell = a is independent of u, and we demanded that 
the fluctuation be fixed at the endpoint in the bulk of AdS space so only the u = e endpoint 
of the u integral contributes. We thus have 



(C ) = a. 



(A.23) 



Arguments very similar to those for eft and {O^) apply also for A z and (J z ). Recalling that 
we absorbed a factor of 2ira' into the vector field above eq. ( 3.10|) , we have 



( J z ) = (2-ira) lim 



1 1 



X osub 
0D D7 



o e q 



The variation of the action is 



r nsub 

db D7 



6d u A z 



X osub 
0D D7 

Sd u A z 



du d u 5A 2 



-P5A z (u 



(A.24) 



(A.25) 



where in the second equality we used the fact that on-shell Sd D J is independent of u and 
in the final equality we demanded that the fluctuation be fixed at the endpoint in the 
bulk of AdS space. Notice that when evaluated on a solution S q d J must be identical 
to the variational derivative of the Legendre transform of the action defined in eq ( 3.14| ): 



5S a j^ 7 b 



Sd u A z - 8d u A 
ingredients, we find 



. The latter is the integration constant j3 in eq. ( [3.16 ). Assembling these 



(J z ) = -(2ttc/)/3. 



(A.26) 



B. Comments on Stability 

In this appendix we do a partial analysis of fluctuations of spinning D7-branes to look 
for indications of instabilities. Since Minkowski embeddings with worldvolume horizons 
are not regular, we focus only on black hole embeddings. We will not compute the full 
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spectrum of quasi- normal modes: we will consider only zero-momentum fluctuations, and 
study only their behavior near the AdS-Schwarzschild horizon. 

To our knowledge the only holographic calculation of quasi-normal modes for spinning 
probe D-branes appears in ref. [69], for a different system. In ref. [69] the field theory 
was M = 4 SYM on a spatial three-sphere 5 3 at zero temperature with massless M = 2 
hypermultiplets and zero magnetic field but finite axial chemical potential. Instabilities 
did indeed appear there which most likely survive the infinite- volume limit in which the S 3 
decompactifies. Our analysis is complementary to that of ref. [69], in the decompactified 
limit, in that we work with nonzero mass, temperature and magnetic field. 

We will use the radial coordinate p defined in eq. ( p.2[ ) and our worldvolume scalar 
dual to O m will be 6(p) as in appendix |A|, which is related to r and R(r) as in eq. ( |A.2j ). 
The boundary is at p — > oo and the worldvolume horizon is at /?* defined in eq. (CT) . The 
equation of motion for 9{p) comes from eq. ( [3.20 ), with r — > p and R(r) — > 9{p). 



To study small fluctuations we will use the full nonlinear equations of motion for the 
embedding that we derived in section || and expand around a background solution to linear 
order. For illustrative purposes we will begin with the trivial background solution, 9{p) = 0, 
which describes massless flavor fields. 

Close to the horizon p* = 1, linearized fluctuations 89(t,r) = e~ lk ° tr d{p) obey an 
equation of motion independent of the value of the magnetic field, 

We can transform this into a Schrodinger form by defining $(p) = (p — l)~ 1 / 2 i/j(p), 

-i>" - J^r§^ = fc V (b.2) 

For any u ^ 0, this equation allows for an infinite set of negative energy /cq < solutions 
representing bound states. This is the same situation as for a field in AdS with a mass 
below the Breitenlohner-Freedman bound, which we take as strong evidence of an instability 
of fluctuations in 9, based on previous experience in similar holographic models [69,97]. 
However, a full analysis of the equations of motion is required to demonstrate the existence 
of an unstable mode. 

We can do a similar analysis for nontrivial black hole embeddings 9{p) ^ 0, although 
the expressions are much more involved. Setting A;o = and expanding close to the 
worldvolume horizon p = p* , we can again derive a Schrodinger-type equation, 

M — - 

{p-p*r 

Here ip is related to i? by a factor that is a power of (p—p*), where the power is a complicated 
function of the parameters to, B, p*, etc. Similarly M is a complicated function of the 
parameters. The condition to avoid negative-energy solutions to the Schrodinger equation 
is M > 0. We plot M as a function of p* for various values of B and several values of 
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(a.) 



(b.) 




(c.) 



Figure 4: The value of M in eq. (B.3) as a function of p 
horizon, (a.) uj = 0.1, (b.) uj = 1, (c.) u 



(d.) 



the position of the D7-brane worldvolume 
2, and (d.) uj = 10, in units of the AdS radius. The red 



(solid), orange (dot-dashed), green (large dashed), cyan (medium dashed) and blue (dotted) curves 
correspond to -B/7 2 = 0.005,0.5,5,50 and 500, where the temperature of the AdS-Schwarzschild 
black hole is T — 7/71". A value M > is highly suggestive that the D7-brane embedding is stable, 
while a value M < is strong evidence that the embedding is unstable. 



positive uj in fig. |] and several values of negative uj in fig. pL With our conventions physical 
embeddings are those with uj > 0, where the body of the D7-brane drags behind the part 
at the boundary. For uj > 0, we find embeddings that satisfy M > 0, for sufficiently small 
uj and B, while for larger values of uj, even solutions close to the trivial one (which has 
p* = 1) have M < 0. As the magnetic field increases, the region where M > shrinks, and 
can disappear completely if the magnetic field is large enough. 

Notice that the behavior of fluctuations does not connect smoothly with the chirally- 
symmetric case = 1. We can see non- analytic behavior at p* = 1 by comparing the 
uj > curves in Figure § with the uj < curves in Figure P the limiting value of M at 
p* = 1 is different for positive and negative uj, despite the fact that the equation of motion 



for fluctuations about the trivial solution, eq. (BJ2), is independent of the sign of uj. 



We should stress again that we cannot decisively conclude whether embeddings are 
stable or not without a full analysis of the fluctuations. We should also mention that we 
have not considered other possible sources of instabilities that can appear for more generic 
fluctuations, for example fluctuations with nonzero momentum. 
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(a.) 



(b.) 



Figure 5: The value of M in eq. (B.3) as a function of p*, the position of the D7-brane worldvolume 
horizon, (a.) uj — —0.1, (b.) uj = — 1, in units of the AdS radius. The red (solid), orange 
(dot-dashed), green (large dashed), cyan (medium dashed) and blue (dotted) curves correspond to 
-B/7 2 = 0.005,0.5,5,50 and 500, where the temperature of the AdS-Schwarzschild black hole is 
T = 7/71- . Notice that the values of M at p* = 1 are very different from those in fig. |J, indicating 
some non-analyticity there. 



C. Effective temperature, entropy and heat capacity of spinning branes 



As mentioned in section 3.1, we can associate a temperature with the worldvolume horizon 
of spinning D7-branes. Recall the induced metric of eq. ( 3.28p , 

ds 2 m = 9udt 2 + 2gg 7 dtdr + gfjdr 2 + g xx dx 2 + g ss ds 2 s s. (C.l) 

Since the induced metric depends only on the radial coordinate r, we can diagonalize it 
with the coordinate transformation 

D7 

di = dt + %~dr, (C.2) 
9tt 



which gives us 
where 



ds 2 D7 = g§ 7 dt 2 + gP/dr 2 + g xx dx 2 + gssds% , (C.3) 



n D7 _ D7 D7 _ D7 _ (fffr 7 ) 2 (n A) 

% — 9tt 1 9ff — 9rr £,1 ■ y^-^l 

9tt 



If g5 7 vanishes at some point r = r*, then we can define a Hawking temperature T# by 
demanding regularity of the Euclidean geometry. Defining the Euclidean time r = it, we 
have 



1 \g p ~ 7 \' 

T ~ T + — — , 1 H ~ 



(C.5) 



Tjj is the Unruh temperature that an accelerating observer on the brane feels. In the 
field theory, the flavor degrees of freedom effectively feel the temperature Tjj while the 
adjoint degrees of freedom remain at zero temperature. We expect two-point functions of 
flavor operators to take a form characteristic of thermal diffusion with temperature Tjj, 
even if the actual temperature is zero [51]. The fact that Tjj can differ from the actual 
temperature is another indication that the system is not in thermal equilibrium. 
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As mentioned in section ^, in the flavor sector only the fermions carry U(1)a charge. 
The scalars are neutral. Our D7-brane solutions with a worldvolume horizon have nonzero 
angular momentum, corresponding to states with a nonzero U(1)a charge, i.e. states 
that naively should have a density of strongly-coupled fermions. An obvious question is 
whether we can detect a Fermi surface. The most direct way to do so would be to com- 
pute two-point functions of fermionic operators, looking for a peak at some finite (Fermi) 
momentum [98-101]. Another approach is to look for the de Haas- van Alphen effect, in 
which the magnetization, magnetic susceptibility, conductivity, and other observables ex- 
hibit oscillations as a function of the magnetic field whose period is inversely proportional 
to the cross-sectional area of the Fermi surface. Detecting such oscillations may require 
computing l/N c effects, as in ref. [102]. Both of these approaches are beyond the scope of 
this paper. 28 

We will follow ref. [106] and characterize our system by asking how the heat capacity 
cy scales with temperature T at low temperatures: a Fermi liquid will have cy oc T while 
a Bose liquid will have cy oc T d in d spatial dimensions. One result of ref. [106] was that 
for D7-branes describing a U(l)y density at low temperatures the heat capacity scaled as 
cy ~ T 6 , which looks like neither a Fermi nor Bose liquid, and hence may be a whole new 
kind of quantum liquid. Here we are asking what happens with a U(1)a density. 

The problem in our case, of course, is that our system is not in equilibrium, so how 
can we define a heat capacity, or any thermodynamic quantity? If the degrees of freedom 
on the D7-brane were gravitational, i.e. if they included spin two fields rather than just 
spin zero and one, then we would have a natural definition of entropy as the area of the 
worldvolume horizon, which would provide a starting point for a thermodynamic analysis. 
Lacking that, we must look for a substitute. The obvious candidate is the on-shell D7-brane 
action, which in equilibrium indeed represents the flavor contribution to the free energy. In 
our non-equilibrium setting, we will define an "effective free energy" as minus the on-shell 
D7-brane action, and extract an entropy and a heat capacity by taking derivatives with 
respect to the effective temperature. Whether these quantities satisfy thermodynamic laws 
is not at all clear. In particular, notice that to, which is equivalent to the axial chemical 
potential, and the temperature Tjj are not independent variables for a fixed \m\. 

The Lorentzian D7-brane action (not just an action density) is 

S D7 = -N J dtdr d 3 x ^-g? 7 gPJg xx (l + ^P) + yBw J dtdr d 3 x A' z {r). (C.6) 

28 Our system, and more generally any system with a CME, has an interpretation in condensed matter 
physics not only in terms of Fermi physics but also in terms of topological insulators. An action of the form 
a(t, x)F A F has been proposed not only as an effective action for the CME [103] but also as an effective 
action for (3+l)-dimensional topological insulators [104]. While a(t,x) = fist produces a CME, an a(t,x) 
independent of t but varying in space as a step-function from zero to tv produces a T-invariant topological 
insulator. An early derivation of the CME in the condensed matter context appears in ref. [1], and an 
analysis of the CME in the topological insulator context appears in section IV B of ref. [104] . Indeed, if the 
QGP created in heavy-ion collisions does in fact exhibit a CME, then the QGP may be a kind of topological 
insulator. A holographic analysis of T-invariant topological insulators using D7-branes in AdS^ x S appears 
in ref. [105]. 
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Here and below we do not write the counterterms of appendix [A], although we include them 
in our calculations in order to obtain a finite, renormalized on-shell action. We can obtain 
the Euclidean action Sp 7 by taking 



S 



D7 



-iS D1 (i 



-IT I 



-s D1 {i 



(C.7) 



Since the D7-brane Lagrangian is time-independent in our case, the replacement t — > r 
acts only on the measure of the integral. Accordingly, the equation of motion in Euclidean 
signature is the same as in Lorentzian signature. Notice that we do not Wick rotate to, so 
the worldvolume scalar c/> becomes imaginary. From a field theory point of view we could 
equivalently begin in Euclidean signature with an imaginary axial chemical potential. 29 

We define the effective free energy density as / = TnSjj 7 /vo1r3. The effective entropy 
density s and the heat capacity cy are then 



df 



dT, 



11 



B,c 



cy = Th 



dT, 



l-i 



B,c 



(C.8) 



In fig. ^ we plot s and cy as functions of Th for cq = 10~ 3 . We find that cy oc Th, similar 
to a Fermi liquid. We also find that s does not vanish in the Th —> limit. Such residual 
zero-temperature entropy, i.e. a large degeneracy of states, is another signal that our 
solutions may not be stable, since a generically any perturbation will break the degeneracy 
and the system will settle in a new, presumably non-degenerate ground state. Notice also 
that our entropy density s is not proportional to the area of the induced horizon, which 
goes as uj^Ol ~ T H for small Th and hence vanishes at Th = 0. 




0.001 0.002 0.003 0.004 0.005 



a. 




0.001 0.002 0.003 0.004 0.005 



(b.) 



Figure 6: (a.) The effective entropy density s versus the effective temperature Th, both in units 
of the magnetic field B, for Co = 10~ 3 , where the flavor mass is \m\ = VX^- in units of the AdS^ 
radius, (b.) The heat capacity cy versus Th, both in units of J5, for the same value of cq as in (a.). 

Calculations similar to the above were done for spinning probe D3-branes without 
a magnetic field in ref. [51]. In that case the on-shell action was not renormalized, and 
divergent terms appeared in the free energy and entropy density proportional to powers 



29 For a probe brane system where an imaginary worldvolume electric field was used to guarantee that 
the Euclidean equations of motion were identical to the Lorentzian ones see ref. [107]. 
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of a UV cutoff. The authors of ref. [51] proposed that the entropy be regarded as an 
entanglement entropy between the flavor sector and the M = 4 SYM sector, in which case 
the UV divergences are the usual ones, with the leading divergence having a coefficient 
proportional to the area of boundary of the entangled region [108, 109] . The UV divergences 
can also be understood from the holographic prescription for entanglement entropy [110- 
112]. We renormalize the action using w-dependent counterterms, so from the perspective 
of ref. [51], we are extracting some part of the finite contributions to the entanglement 
entropy. 
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